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Chapter 1 


Linear Partial Differential 
Equations 


Problem 1. Show that the fundamental solution of the drift diffusion 
equation 


du Pu _ OU 
dt Ax? | Ox 
is given by 
_ 1 (x — 29 + 2t)? 
u(a,t) = exp ( Tia a 


Solution 1. 


Problem 2. (i) Show that 


pnpy= ot () 

(ii) Show that 
Dy (f+ 9) = (-1)" Dz (9: f)- (2) 

(iii) Show that 
D(f-f)=0, for m odd (3) 


Solution 2. 


2 Problems and Solutions 


Problem 3. (i) Show that 
D7? Di (exp(kix — wit) - exp(kex — wet)) = 
(ky — ko)” (—w1 + we)” exp((ki + ko) — (wi + we)t) (1) 
This property is very useful in the calculation of soliton solutions. 
(ii) Let P(D;, Dz) be a polynomial in D; and D,. Show that 
P(Dz, Dt) (exp(kix — wit) - exp(kex — wet)) = 


P(ki — kz, —w1 + w2) 
P(ky + ko, =O we) 


P(Dz, Di) (exp((ki + ka)u — (wi +w2)t)- 1) (2) 


Solution 3. 


Problem 4. Consider a free particle in two dimensions confined by the 
boundary 


G:= { (x,y) : |ey| = 1}. 


Solve the eigenvalue problem 


Ay + k*y =0 
where 
Gone 2nE 
h2 
with 
Ya =0 
Solution 4. 


Problem 5. Consider an electron of mass m confined to the x — y plane 
and a constant magnetic flux density B parallel to the z-axis, i.e. 


B=] 0 
B 


The Hamilton operator for this two-dimensional electron is given by 


T= (p + eA)? = 1 


a 2 a 2 


where A is the vector potential with 


B=VxA 
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and 


or 
—yB 
A= 0 
0 
(ii) Use the second choice for A to find the Hamilton operator A. 
(iii) Show that ; 
[H, pe] = 0. 


(iv) Let k = p,/h. Make the ansatz for the wave function 
o(x,y) = e**o(y) 


and show that the eigenvalue equation H yw = Ew reduces to 


( a | a (y w)) o(y) = Edly) 


where 
_ eB __ hk 
mane, EB 


(v) Show that the eigenvalues are given by 


1 
Bn = (m4 5 ) hie, n=0,1,2,.... 


Solution 5. (i) Since 


OA, /Oy — OA, /Oz 
VxA:= | 0A,/02 — 0A,/0x 
OA, /Oxz — OA, /Oy 


we obtain the desired result. 
(ii) Inserting A, = —yB, Ay = 0, Az = 0 into the Hamilton operator we 
find 


Dr B)?. 
_, (Px — eyB) 


4 Problems and Solutions 


(iii) Since coordinate x does not appear in the Hamilton operator H and 


it follows that [H,p,] = 0. 
(iv) From (iii) we have 


pev(a.y) = ih i(e,y) = Akw(a,y), — AW(2,y) = EY(2.y). 


Inserting the ansatz ~(x,y) = e'**w(y) into the first equation we find 
b(x,y) = e™*9(y). 
(v) We have 
na 1 ~ ~ \ ika 
y= In | (Pe — eyB)” + py)e™* $(y) 
85 F 
= 5 e'*((hk — eyB)” + By)o(y) 
m 
It é 
= ae (eB — hk/eB)” + p5)0(y) 


‘, 2 Ae 
= eit (“2 (y— yo)? + ji oly). 


Now the right-hand side must be equal to Ey)(x,y) = Ee’**(y). Thus 
since 


the second order ordinary differential equation follows 


( ois a (y w)) o(y) = Egy). 


(vi) The eigenvalue problem of (v) is essentially the harmonic oscillator 
except the term is (y — yo)? instead of y?. This means that the centre of 
oscillation is at y = yo instead of 0. This has no influence on the eigenvalues 
which are the same as for the harmonic oscillator,namely 


1 
Ey = (145) te, n=0,1,2,.... 


Problem 6. Consider the Schrodinger equation 
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Find the coupled system of partial differential equations for 


p=, HS) (=) ; 


Solution 6. First we calculate 


6) 7) Ow* 
Fe ie) = et 
Inserting the Schrodinger equation and 
noe _ 1 2 
yields ; 
Pp 1 
OE. - 5 We) dima oY , 


We set w = e®)+#S(® where R and S are real-valued. Next we separate 
real and imaginary part after inserting this ansatz into the Schrodinger 
equation. Then we set v = VS. 


Problem 7. Consider the conservation law 


Oc(a, t) ps3 Oj (x,t) 
Ot Ox 


= 0 


where Ae(-r,t) 
Lr, 

5(z,t) = —D(x) 

x 


The diffusion coefficient D(a) depends as follows on x 


D(e)=Do(l+gx) 9 0. 


Thus dD(x)/dx = Dog. Inserting the current j into the convervation law 
we obtain a drift diffusion equation 


Oc(a,t) Oc(a, t) 0 c(a, t) 
Bp ag. OE eg 


The initial condition for this partial differential equation is 
c(a,t = 0) = Mod(x) 
where 6 denotes the delta function. The boundary conditions are 


j(a@ = 2, t) = j(a = +00, t) = 0. 


6 Problems and Solutions 


Solve the one-dimensional drift-diffusion partial differential equation for 
these initial and boundary conditions using a product ansatz c(#,t) = 
T(t)X (a). 


Solution 7. (Martin) Inserting the product ansatz into the one-dimensional 
drift diffusion equation yields 


1 dT(t) 1 dX(x) 1 d2X(x) 


Te a SRG as 


Setting the left-hand and right-hand side to the constant — Dog? /4, where 
x is the the speration constant (x > 0) we obtain for T(t) the solution 


T(t) = C(x) exp (--41) 


where the C(x) has to be calculated by the normalization. For X(x) we 

find a second order linear differential equation 

d?X dx 2 
XI xX 


dx |! dx 4 mae 


(1 + gx) 


Introducing the new variable u = 1+ gx and Y(u) = u-!X(u) we obtain 


wy dY x 
2 { — 
U du2 3u om (a+ *u) Y =0. 


The solution of this second order ordinary differential equation is 


¥(u)=— Jove, X(u) = u¥(w). 


Here Jo(,/xu) is the Bessel function of the first kind of order n = 0. We 


have a p 
In() = » oan a 
Thus re , ms 
Jota) = > ar () 


It follows that 
X(u)T(t) = Clo.) Jol Fa) exp (--41) | 


Now we integrate over all values of the separation constant y > 0 


sas) = [aco RcEae exp (-22er) 
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The normalization constant C(x) is determined as follows. Inserting t = 0 
(initial condition) we obtain 


det =0) = f° dyC(x)so Vx 0) 


Multiplying with [ dxJo(\/x/(1+ gx)) and using the orthogonality of the 
Bessel functions we obtain 


1 
C(x) = qMogJo( Vx). 
It follows that 


cast) = “PE [ay Jol yX)J0l VC Fe) exr(—xDoa?t/4) 


The y integration can be performed and we obtain 
th= if /14 pa 
c(z, t) Dogt ( Dogtt a) exp ( Dow? *) 


where Jo is a modified Bessel function of order 0. J, is defined by 


a oy 1 ay U+2k 
Leet Tia) = > are aD) 3) 


For g — 0 we obtain 


Goa 
Cx — x é 
ye ia \. ADyb 


Problem 8. Consider the time-dependent Schrédinger equation 


_O 1 
igv=(-3A+V@)) 0 (1) 
where . . ‘ 
0 0 0 
A= at zt 2 
Oxi 0x5) ONS 
Let 


8 Problems and Solutions 


Show that the time-dependent Schrédinger equation can be written as the 
system of partial differential equations (Madelung equations) 


Op 


ap ewe) = (2 ee ase) ) (2) 
oe 4(v-V)jv=-V (vin) = ae) (3) 
Solution 8. To find (2) we start from (1) and 
ne = ( a | v(x) wv. (4) 
Now from p = ww* we obtain 
a BEY 


Inserting the time-dependent Schrédinger equations (1) and (4) we obtain 


ee (-a + v(x) Ww (-a + v(x) yr 


at ih 2m 
i ae en Py ae . 
= ih (- awe ) GE (-vaw ) 
al 1 . 3 Ow Ov* 
=a \~a, | 40 +8) Be, 8a; 


Equation (3) can be derived by writing w in to so-called Madelung form 
(x, t) = exp(R(x, t) + iS(x, t)) 

where R and S are real-valued functions. Substituting this expression into 

the Schrédinger equation, dividing by w, separating real and imaginary 


part, and taking the gradient of the equation in S, where v = VS. The 
Madelung equations are not defined on the nodal set since we divide by w. 


Problem 9. Consider the time-dependent Schrédinger equation 


Ov) 
ih a —5— Vw(x, t) + V(x, t)a(x, t). 


Consider the ansatz 


w(x, t) = (x, t) exp(tmS(x, t)/h) 
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where the functions ¢ and S are real. Find the partial differential equations 
are @ and S. 


Solution 9. Since 
Oy _ 
Ot 
and 
Ow = 


= 
Ox; 


we obtain the coupled system of partial differential equations 


0 2 2 = 
ae + V+ (VS) =0 


a 1 (OC -(R?/2m)V2¢ 
HVS t(VS-V)VS= (v ‘ 


+ vw) ; 
m 
This is the Madelung representation of the Schrodinger equation. The term 


ep eLUNe 
b 


of the right-hand side of the last equation is known as the Bohm potential 
in the theory of hidden variables. 


Problem 10. Consider the Schrodinger equation HW =EWV ofa particle 
on the torus. A torus surface can be parametrized by the azimuthal angle 
¢ and its polar angle @ 


x(@, 0) = (R + acos(9)) cos(¢) 
y(%, 0) = (R + acos(#)) sin(¢) 
2(@, 6) = asin(0) 


where R and a are the outer and inner radius of the torus, respectively such 
that the ratio a/R lies between zero and one. 

(i) Find H. 

(ii) Apply the separation ansatz 


V(0, d) = exp(img)y(4) 


where m is an integer. 


Solution 10. (i) We find 


Fi hs oP sin(6) an 1 oe? 
~ 2 \a?2 00? a(R+acos(@)) 00 ' (R+acos(6))? 0¢? 
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(ii) Let a = a/R. We obtain the ordinary differential equation 


ay asin(0) diy ma? 
dé? 1+ acos(@) d@ = (1+ acos(6))? 


¥(8) + By =0 
where 6 € [0,27], a € (0,1), and @ is a real number with 6 = 2Ea?. 


Problem 11. The linear one-dimensional diffusion equation is given by 


du _ you 

Ot Ox? 
where u(x,t) denotes the concentration at time ¢ and position x € R. D 
is the diffusion constant which is assumed to be independent of x and t. 
Given the initial condition c(z,0) = f(x), « € R the solution of the one- 
dimensional diffusion equation is given by 


t>0, -w<4<o 


u(x,t) = he. G(a, ta’, 0) f(x’)da’ 


where 


T ogl\y __ 1 (x = a!) 
rn) = meme ( 4D(t — ) . 


Here G(x, t|z’, t’) is called the fundamental solution of the diffusion equation 
obtained for the initial data 6(a — x’) at t = t’, where 6 denotes the Dirac 
delta function. 

(i) Let u(x, 0) = f(x) = exp(—2x?/(2c)). Find u(z, t). 

(ii) Let u(v,0) = f(a) = exp(—|z|/o). Find u(z, t). 


Solution 11. 


Problem 12. Let f : R®? — R® be a differentiable function. Consider the 
first order vector partial differential equation 


Vxf=hkf 


where & is a positive constant. 
(i) Find V x (V x f). 
(ii) Show that Vf = 0. 


Solution 12. 
Problem 13. Consider Mazwell’s equation 


ae Yee. divE=0, divB =0 


Wie Ba OE! Ot 
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with B = woH. 

(i) Assume that FE. = E3 = 0 and By = B3 = 0. Simplify Maxwell’s 
equations. 

(ii) Now assume that FE, and Bz only depends on «3 and ¢ with 


Ey (x3, t) => f(t) sin(k3x3), Bo(x3, t) = g(t) cos(k32'3) 


where kg is the third component of the wave vector k. Find the system of 
ordinary differential equations for f(t) and g(t) and solve it and thus find 
the dispersion relation. Note that 


OE3/O0x2 = OE2/0x3 
VxE= OF; /0x3 — OE3/0x1 
OE2/O0x, oe? OE; /0x2 


Solution 13. (i) Maxwell’s equation written down in components take 
the form 


OB, OE: OFs 


ot 0x3 Ore = 

OBo OF3 OF, 9 

Ot Or, 3023 

OB; OE, OE» 

Ot 0x2 = Oxy =e 
1dE, OB; | 2Ba _, 
c2 Ot Oro | Ox3 
1B, OB | OBs _, 
c2 Ot Or3 Oz, 
1 dE; OB2 | Bi _, 
c2 Ot Or, Oxy 

and 
OF, OF» OF3 OB, OBo OB; 


| 
T 


='0. 


=0 
0x4 0x2 i 0x3 ‘i Oxy ie 0x2 = 0x3 
With the simplification FE. = E3 = 0 and B, = Bs = 0 we arrive at 
OE OB OE OB OE OB 
—t=0, —24+—*=0, —*=0, +2 =0 
0x2 Ot 0x3 Ox ot 0x3 
and 


OE,  OBy 

Oxy , Ox 
(ii) With the assumption that E, and Bz only depend on x3 and ¢ the 
equations reduce to two equations 


OB. OB, OB, . 0B, _ 
Oe be SO ba 
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Inserting the ansatz for E, and Bo yields the system of differential equations 
for f and g 
dg 
2 
C R39 ’ dt a 3h 


These are the equations for the harmonic oscillator. Starting from the 
ansatz for the solution 


f(t) = Asin(wt) + Bcos(wt) 


where w is the frequency, yields 


Aw Bu . 
g(t) = Oks cos(wt) — Dh, sin(wt) 
and the dispersion relation is 


a eke: 


Problem 14. Let k be a constant. Show that the vector partial differen- 
tial equation 
Vxu=ku 


has the general solution 
1 
u(21, 22,03) = V x (cv) + ra x V(cv) 


where c is a constant vector and v satisfies the partial differential equation 


V-u + ku = 0. 


Solution 14. 


Problem 15. Consider the partial differential equation of first order 


() 3) ) 
(2 Ox a) io 


Show that f is of the form 


{Gx Pasta — Bye a): 


Solution 15. We set Yi2 = T%1 — &Q, Y23 = T2 — V3, Y31 = V3 — X34. Then 
applying the chain rule we have 


Of Of Oyi2 Of Oys1 Of Of 


dx, Oyi2 Ox, © Oygi Oxy ~— Oya ~— Oysi 
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Analogously we have 


af 9 a af a a 
dry O 0 Ox, 0 O 


Thus f(x — x, — a,x — x) satisfies the partial differential equation. 


Problem 16. Consider the operators 


0 (a) 0 
x + a; => ,t)5-. 
Det ppt Dem ae) a. 
Find the commutator 
[Dz, Di) f(a, t) 
Solution 16. We obtain 
Ou 
[Dz, Dilf = age 
v 


Thus 
[Dz, Di) _ UgDre. 


Problem 17. Let C be a constant column vector in R” and x be a 
column vector in R”. Show that 


Cex C= Vero 


where V denotes the gradient. 


Solution 17. Since 


x™C = >> 2,C; 
j=l 
and 3 
oe T = 
ee C=C; 


we find the identity. 


Problem 18. Consider the partial differential equation (Laplace equa- 
tion) 
Ou Ou 
da? * Oy? 
with the boundary conditions 


=0 on [0,1] x [0,1] 


u(z,0)=1, u(#,1)=2, u(0,y)=1, u(l,y) =2. 


14 Problems and Solutions 


Apply the central difference scheme 


(53) re Uj—1,k — 2U5,k + Uj+1,k (=)  Uik-1 a QU; k + Uj,.k+1 
Ox? ik (Az)? : Oy? pk (Ay)? 


and then solve the linear equation. Consider the cases Ax = Ay = 1/3 and 
Aa = Ay = 1/4. 


Solution 18. Case Ar = Ay = 1/3: We have the 16 grid points (j,k) 
with j,k = 0,1,2,3 and four interior grid points (1,1), (2,1), (1,2), (2,2). 
For 7 = k = 1 we have 


Uo,1 + U21 + U1,0 + U1,2 — 4u1,1 = 0. 


For 7 = 1, k = 2 we have 
Uo,2 + U2,2 + U1,1 + U1,3 — 4u1,2 = 0. 
For j7 = 2, k = 1 we have 


U11 + 3,1 + 2,9 + U2,3 — 4U2,1 = 0. 


For 7 = 2, k = 2 we have 


U1,2 + U3,2 + U2,1 + U2,3 — 4U2,2 = 0. 


From the boundary conditions we find that 
Uo,1 = 1, U1,0 = 1, U0,2 = 1, U1,3 = 2 
Uz1i=2, Uwo=l, uwz2=2, Uu23 = 2. 


Thus we find the linear equation in matrix form Au = b 


—4 1 1 0 U1,1 —2 
1 —4 0 1 U1,2 = —3 
1 0 —4 i} U2,1 — —3 
0 1 1 —4 U2,2 —4 


The matrix A is invertible with 
7/2 1 1 1/2 
mes 1 7/2 1/2 1 
12 Le Dp 2e FAQ" 
1/2 1 1 7/2 


Al = 


and thus we obtain the solution 


U1,1 5/4 
_ U1,2 _ 3/2 
aE U2,1 _ 3/2 


U2,2 7/4 
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Case Ax = Ay = 1/4. We have 25 grid points and 9 interior grid points. 
The equations for these nine grid points are as follows 


Uo,1 + U2.1 + U1,0 7 U1,2 — 4uy4 =0 
Uo,2 + U2,2 + U1,1 + U1,3 — 4U1,2 =0 
Uo,3 + Ua2,3 + U1,2 + U1,4 — 4t1,3 =0 
U1,1 + U3,1 + U2,9 + U2,2 — 4u21 =0 
U1,2 1 U3,2 T U21 T U2,3 — Aus 2 =0 
U1,3 + U3,3 + U2,2 + U2,4 — 4U2,3 =0 
U2,1 + U4, + U3,0 + U3,2 — 4u3,1 =0 


U2.2 + U4,2 + U3,1 + U3,3 — 4u3,.2 = 0 


U2,3 + Ua4,3 + U3,2 + U3,4 — 4U3,3 = 0. 
From the boundary conditions we find the 12 values 
ui =1, uo =1, vo2=1, uo3=1, wi4=2, U20=1, 


U24=2, U41 =2, Ugo =1, Usa = 2, Uag = 2, Uz = 2. 


Inserting the boundary values the linear equation in matrix form Au = b 
is given by 


—4 1 0 1 0 0 0 0 0 ut, —2 
1 -4 1 0 1 0 0 0 0 U1,2 = 
0 1 -4 0 0 1 0 0 0 U1,3 —3 
1 0 0 -4 1 0 1 0 0 U2, —1 
0 1 0 1 —4 1 0 1 0 U2,2 = 0 
0 0 1 0 1 -4 0O 0 1 u2,3 2. 
0 0 0 1 0 0 -4 1 0 Us, —3 
Oo AG 0 Gr: > A: a tee “a U3, = 
0 0 0 0 0 1 0 1 —-4 U3,3 —4 


The matrix is invertible and thus the solution is u = A~!b. 


Problem 19. Consider the partial differential equation 


OF of 5 
(PB) yarpno 
(i) Apply the transformation 


2 


u(x, y) SS ie v(x, y) = 2ry, W(x, y) an w(u(x,y), v(@,y)): 


(ii) Then introduce polar coordinates u = rcos¢, v = rsin@. 


16 Problems and Solutions 


Solution 19. (i) Note that 
u? + uy? = (a? + y7)?, ety = /y2 + y2, 
Applying the chain rule we have 


Op _ dpdu  dpdv _, ab 
dx OuOx | «~OVOX "Ou Foy 


and 
Py db | 20y Py 20y 
Ox? * os oe Ou? Pee Nan + ay Ov 
Analogously we have 
Op _ dpdu  dbdv__, Ob 4, OW 
Oy Oudy  OdOvdy — ae Ov 
and % ts a zs 
Py do aH ab Od 
Oy2 Fre Ou? a ee Ov2- 


Consequently we find the partial differential equation 


Introducing polar coordinates we have 
af ob\ 10% PK 
Tr ¥ oy + i = 0. 
Or \ Or r 0g? 4 


One can now try a separation ansatz 


vr, @) = F(r)g(%) 


nae ae, 1 Pg 
Ff) dr (-f) “edge a 


Thus we find the two ordinary differential equations 


r (df ae ee 
+(f +s) wie Bese 


which yields 


1d’g 
g dg? 


where C is a constant. 
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Problem 20. Starting from Maxwell’s equations in vacuum show that 


kx E=woB. 


Solution 20. 


Problem 21. Show that the linear partial differential equation 
Pu Pu 
06? Ot? 
admits the solution u(t,@) = sin(@) sin(t). Does this solution satisfies the 
t,0 0 


boundary condition u(t, = 7) = 0 and the initial condition u(t = 0,0) = 
0? 


Solution 21. 


Problem 22. Consider the Hamilton operator 


Beth Be ae 
H= d)? = 90)" 
5 Din, ag a 5 (a2 an ) + (93 — @ I) 


where d is the distance between two adjacent atoms. Apply the linear 
transformation 


€(q; 92,93) = (42-4) —d 
n(q1, 92; 93) = (93 = q2) —d 
1 


X(41, 4%; 93) = (M141 + Moq2 + M3q3) 


M 


where M := m,+m2+mz3 and show that the Hamilton operator decouples. 


Solution 22. The Hamilton operator H takes the form H = Hom + An, 
where 


: fy 
Hom = aye pu := —tha > 


* h? 0? he io? ho k 
Ay = ; 5+ + ~(€7 + n°) 
2p OF? = 2a. ON? M2 OEON 2 
where py = (mmz2)/(m1 + m2), 2 = (Memz3)/(m2+m3). Thus the centre 
of mass motion is decoupled. 


Problem 23. Consider the Hamilton operator for three particles 


* 1 0? o? oO? 
(soz * ang + dg) ~ 8060 — 220822 —29) 


2m 
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and the eigenvalue problem Hu(21, x2, 3) = Fu(a1,%2,x73). Apply the 


transformation 
2/1 
Hi (21; 29; 23) = 3 5 ("1 +22) — v3 


1 
Y2(@1, 2, £3) = va = £2) 


(x1 +%2+ 43) 


Co] rR 


y3(X1, £2, 13) oa! 
U(y1 (1, 2, £3), Yo(%1, T2, €3), Y3(X1, 2, %3)) = U(#1, £2, V3) 


where y3 is the centre-of mass position of the three particles and y, y2 give 
their relative positions up to constant factors. Find the Hamilton operator 
for the new coordinates. 


Solution 23. First we note that 


O(a@1, La, £3) _ 


O(y3, V1, Y2) 


= 


The Hamilton operator takes the form 


~ 1? 1 f/f oP 


Consequently the variable can be ys can be separated out. 


Problem 24. Consider the four dimensional Laplace equation 


Oru 3 Pu a Pu Pu 
Oxi One OES Y Ba 


Transform the equation into polar coordinates (7,0, ¢, x) 


x1(r, 0, ¢, x) =r sin(@/2) sin((¢@ — x)/2) 
xa(r, 0, 6, x) =r sin(@/2) cos((¢ — x)/2) 
x3(r, 0, 6, x) =7r cos(6/2) sin((¢ + x)/2) 
xa(r, 0, 6, x) =7r cos(0/2) cos((d + x) /2) 

u(x) = u(r(x), A(x), oC), x(x) 


where x = (#1, 2,3, %4). 
Solution 24. In polar coordinates we find 


Vyti + = {V3 ua =0 
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where 
10.20 
Dette eg 
Ve r3 ar’ Or 
0 0 1 fo oP O? 
2 = si 6 
Va sin 0A, sin Op sm? 8 te De 2.cos sx) 


Problem 25. A nonrelativistic particle is described by the Schrédinger 
equation 


~ 2 
Pp = OY 
(= a V(x, 0) w(x, t) =ih at 
where p := —7ihV. Write the wave function w in polar form 


w(x, t) = R(x, t) exp(iS(x, t)/h) 


where R, S are real functions and R(x,t) > 0. Give an interpretation of 
p=. 


Solution 25. The wave equation is equivalent to a set of two real differ- 
ential equations 


2 
7p +V (=) =0 
where ee 
aes oR 


The equation for the time evolution of p is a conservation equation that 
provides the consistency of the interpretation of p as the probabilty density. 


Problem 26. The time-dependent Schrodinger equation for the one- 
dimensional free particle case can be written in either position or momen- 
tum space as 


2 Plant), OU(x,t) 


2m Ox? Ot 


Consider the momentum space approach with the solution 


(p, t) = $o(p) exp(—ip?t/2mh) 
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where ¢(p,0) = ¢o(p) is the initial momentum distribution. We define 
& := ih(0/Op). Find (#)(t). 


Solution 26. We have 
(a(t) = / _# (not 0) a 
=| e(eieo(ryan+ =f rloo(n.t)Pap 
(p)ot 


= (£)o+ a 


Problem 27. Consider the partial differential equation 


2 


0 
Oxdot 


(In(det(I,, + tD£(x))) = 0 


where Df(x) (f : R” — R”) is the Jacobian matrix. Find the solution of 
the initial value problem. 


Solution 27. The general solution is 
det(I, + tDf(x)) = g(x)h(t) 


where g: R° ~ R,h: R—-R. At t = O since det(I,) = 1 we have 
g(x) =h7'(0). Thus 


det(I, + tD£(x)) = ait 


m= 


Thus 


a 
Fy det(In + tDE(x)) = 0. 


Problem 28. The partial differential equation 


du \? Pu \? _ ( Ou 5 Pu \? _ ( ou ; Pu \? 
(7) (saa) (55) (saa) (se) (sas) 
& Ou Ou Ou be Ou Ou Ou O?u _ Ou du Ou Au ) 
Ox, Ore Ox10x3 Ox20r3 | OX, Or3 OX1OX2 OX2OT3 | -OX2 OX3 OX1OX OL OL3 
Ou Ou Ou Ou 


Oxy 0x2 0x3 021022023 i 
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has some link to the Bateman equation. Find the Lie symmetries. There are 
Lie-Backlund symmetries? Is there a Legendre transformation to linearize 
this partial differential equation? 


Solution 28. 


Problem 29. Consider the two-dimensional heat equation 


Ol arr 4 OT 
Ot = Ox? Oy? 


Show that 


Tes yst) = exp ( a) 1 ( e+ uen) 


satisfies this partial differential equation. 


Solution 29. 


Problem 30. Let « > 0. Consider the Fokker-Planck equation 
Ou uO 
€ 


with u > 0 for all z. 
(i) Find steady state solutions, i.e. find solutions of the ordinary differential 
equation 


uO 
s+ =0. 
“Ox? = Ox ei 
(ii) Find time dependent solution of the inital value problem with 


1 
V 27700 


with og > 0 the variance and fio the mean. With o = 0 we have the delta 
function at 0. 


u(t = 0,2) = N,,(z — uo) = e~ (®—Ho)*/ (20) 


Solution 30. (i) We have 
d {du xu 
fae et | 0). 
dx (=F € ) 


dus xu 
ie ee 


Integration yields 
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where C; is a constant of integration. Consequently 


yee? lt22)) = Cer /29), 


dx 


Integrating once more yields 
u(x) = (c | e [28 gg 4 cr) e7® /(2e) 
0 


where C2 is another constant of integration. The condition u > 0 implies 
C = 0 and the condition 


implies 
u(x) = 1-2 /(26), 
QTE 
Thus this time-independent solution is the normal distribution with mean 
0 and variance e. 


(ii) We obtain 


with 
u(t) = poe *, a(t) =€+(o9—€)e7 


Problem 31. Consider the Schrodinger equation 


N92 
— age + 20d, Ha; - te) u(x) = Eu(x) 


j<k 


describing a one-dimensional Bose gas with the d-fnction repulsive interac- 
tion. Show that for N = 2 the eigenvalue problem can be solved with a 
exponential ansatz. 


Solution 31. For N = 2 we have the unnormalized wave function 


u(21, £2) = eikivittkew2 4 e012 ptkawitikixe 


for 21 < x2 where the phase shift 612 is given by 


2 JOS i(ky = kg) 
c+ i(ky —s ka) 


e912 _ 


The wave function for x, > x2 is obtained by using the symmetry relation 
u(@1, 22) = u(w@e,21). For ky = kg we have u(x1, 22) = 0. 
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Problem 32. (i) Show that the power series 
co 
P(t, A) = S- PnX" cos(/nt) 
n=0 
satisfies the linear partial differential equation (diffusion type equation) 


O° PEt) a) 
a2 Az Pt A) 


with boundary conditions 


OP(0, d) 


Fe =O Plts0) = 0(0). 


P(O,A) = DJ pnd” = pl); 
n=0 
(ii) Perform a Laplace transformation 
P(z,X):= / e *' P(t, \)dt 
0 
and show that P(z,A) obeys the differential equation 
Hes as 

2° P(z, X) — zp(A) = —-A—P(z,d) 


On 


with boundary condition 


and 


Poe.) = f exp(—2? In(A/x)) eae. 


Solution 32. 


Problem 33. Consider the Schrédinger equation (eigenvalue equation) 
for an identical three-body harmonic oscillator 


nS eS 
Baus? So tna ae 
2M 4 2 eer 


where M denotes the mass and A,, is the Laplace operator with respect 
to the position vector r;. In the center-of-mass frame the position vector 
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r; is replaced by the Jacobi coordinate vectors R, x, y 


1 
R=-—=(r1+4ret+r3) =0 


v3 
1 23 
x= V6 (2r1 TQ r3) = Vat 


1 
= —(|ro-7r 
y Ta! 2 3) 
V(R(z;), x(rj), y(4j)) = (rr, v2, 43). 
Show that the eigenvalue equation takes the form 


2 


under this transformation. 


Solution 33. 


Problem 34. Show that 


u(a,t) = 


1 
ex 
Vandi? ( ADt 
satisfies the one-dimensional diffusion equation 


Ou D O7u 


Ot ~Ox? 
with the initial condition 


u(x,0) = d(x — 20). 


Solution 34. 


Problem 35. Consider the integral equation 


oats, t) = —ac(a,t) + 2 | c(y, t)dy. 


Let s > 0. Show that 
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with c(a,0) = e~** satisfies the integral equation. Show that taking the 
derivative with respect to x we obtain the linear partial differential equation 


O7c = re 
Otox Ox 


Solution 35. 


Problem 36. Consider the one-dimensional diffusion equation 


Ou’ Ou! 
Ot =) 5 u' (x,t = 0) = u(z) 


Show that 
tt (3t)= — | u(a + 2VDis)e~* ds, x>0 
is a solution. 


Solution 36. 


Chapter 2 


Nonlinear Partial 
Differential Equations 


Problem 1. Consider the system of quasi-linear partial differential equa- 
tions 


Ou Ov | Ov Ou 


OE Oe (Ok Oe 
Ou OV _ 
Oy Ox 


Show that this system arise as compatibility conditions [L,M] = 0 of 
an overdetermined system of linear equations LW = 0, MW = 0, where 
W(2,y,t, is a function, » is a spectral parameter, and the Laz pair is 
given by 


Solution 1. Using the product rule we have 


[L, MW = L(MW) — M(LW) 


9 ONO BD 
=(5 "Be x) Gane az) y 


Ou OW Ou OW Oudv dvdAYV Ov OW Ov OW 


—O.0n Ob n UUs  OuOn | (eos. OR On 
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_ (du Ou oe Ov Ov Ov\ OW 
~ Lat Ox Oy Oy Ox Ox } Ox 


(F Ou Ov Ov ($ =)) ow 
= Vv t + U r t 


Ot Ox Ox Oy Ox” 


Thus from [L, M]V = 0 the system of partial differential equations follows. 
Problem 2. Find the traveling wave solution 
u(x,t) = f(x — ct) c = constant (1) 


of the one-dimensional sine-Gordon equation 


Ou du ; 
Age. = OE = sin(w). (2) 
Solution 2. Let 
s(a,t) =a — ct. (3) 


Then we find the ordinary differential equation for the function f 
d2 
(1- eyet = sin(f). (4) 


Solutions to (4) can be written as elliptic integrals 


f dz 8 
/ \/2(E — cos z) = V1-¢ ©) 


where £ is an arbitrary constant of integration. From (5) it follows that 


f(s) = cos! (2e” ( (6) 


wise)» 


where cd is an elliptic function of modulus y = 2/(E + 1). 


Problem 3. Show that the nonlinear nondispersive part of the Korteweg- 
de Vries equation 


Ou Ou 
= 1 
a + (a+ Bul =0 (1) 
possesses shock wave solutions that are intrinsically implicit 
u(x,t) = f(x — (a + Bu(x, t))t) (2) 


with the initial value problem u(x,t = 0) = f(x), where a, GER. 
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Solution 3. Since 


se = ( a— Bu tse) 


Ou Ou\ ,, 
Se (1-005) f 
where f’ is the derivative of f with respect to the argument we obtain 
Ou\ Ou Ou\ Ou 
(1 at) at ={ ooue att) On (4) 


From (4) it follows that 


which is (1). 


Problem 4. Consider the Korteweg-de Vries-Burgers equation 


UU ns Ou 07 u Pu ey, (1) 
Of op | One One F 


where a1, a2 and a3 are non-zero constants. It contains dispersive, dissipa- 
tive and nonlinear terms. 
(i) Find a solution of the form 


bi 


u(a,t) = (1 + exp(ba(a + bat + b4)))2 


(2) 


where 61, bo, bz and by are constants determined by aj, a2, a3 and a4. 
(ii) Study the case t > co and t > —ow. 


Solution 4. (i) This is a typical problem for computer algebra. Inserting 
the ansatz (2) into (1) yields the conditions 


bs = —(a4b1 + agb2 + a3b>) (3) 
and 


(2a1b1+3a2b2+9a3b3) exp((a1b,b2 t abs t a3b3)t) t (a,b, t 3a2b2 3a3b3) 0. 


(4) 


The quantity b4 is arbitrary. Thus 


2a,b, = —3agb2 — 9a3b3, a,b, = —3agb2 + 3a3b5. (5) 
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Solving (5) with respect to b; and by yields 


12a? a2 
b = — 2 b =o 6 
25a1a3° : 5a3 ( ) 
Inserting (6) into (3) gives 
- 6 a (7) 
3 25 a3 


Thus the solution is 


12a3 ag 6a3 ie 
= 14 t t+b , 8 
CO Ries re ( iad (2 (« hag (8) 


(ii) Two asymptotic values exist which are for 


t — —oo, ul = — (9a) 


and for 
t > oo, urr = 0. (9b) 


Problem 5. The Fisher equation is given by 


Ot Ou? 


where the positive constant a is a measure of intensity of selection. (i) 
Consider the substitution 


v(a,7T(t)) = aula t) T(t) = 5t. (2) 


(i) Show that (1) takes the form 


(ii) Consider the following ansatz 


= S$ > u;(2,t)¢?-?(a, t) (4) 


j=0 


is single-valued about the solution movable singular manifold ¢ = 0. This 
means p is a positive integer, recursion relationships for v; are self-consistent, 
and the ansatz (4) has enough free functions in the sense of the Cauchy- 
Kowalevskia theorem. 

(iii) Try to truncate the expansion (4) with v; = 0 for j > 2. 
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Solution 5. (i) Using the chain rule we find that substituting (2) into (1) 
provides (3). 
(ii) By substituting (4) into (3), we have p = 2 and 


Ovj—2 . Od = 0? v;_2 . Ovj-1 Od ‘ ; Oro 
: (= HG 3) <2) Ox? +2(5—3) Ox Ox AIS Spa Ox? 
ad 2 J 
t © 3) @ 2)u; be + avuj—2 — 6a S- Uj—kUk (6) 
k=0 


for all 7 > 0, where v; = 0 for 7 < 0. For 7 = 0, we have 


wat (2). (7) 


Thus using (7), (6) turns into 


(j- 6) + Lu, (32) a5 (2 t(j 3)u5150 | Puja 


Ox OT Ox? 
j—1 
gn O0j-106 ao . 
2(j — 3) Aa Bi ©) 3)¥j-1 573 + avj-2+ 60 )_¥j-# (8) 
for 7 > 2. The resonance points are 7 = —1 and 7 = 6. The point —1 


corresponds to the arbitrary singular manifold function ¢, while the point 
6 corresponds to the free function vg. For 7 = 1, we obtain from (6) that 


nai (eo), (9) 


~ a\dr Ox? 


We find all functions of v; using (8). 
(iii) Using (6), we let v; = 0, for 7 > 3. Then 


v=szt—4+ 0. (10) 


Hence (8) gives 


mat lll (Gen)? 1 (done) Ger, 1brebe 1 ($2) 
ane ber 4 be 3 dx dx ‘3 oe 12 bx 


(11) 
bea 7a br as 1 Pow Ve ob; 1 1 ear 5 Qrr \ 1 rane 
bs =a bs pea 2 be 2b, 12 de on 
and ‘ 2 
5 = sae + ave(1 — 6v2) (13) 
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where we set ¢, = 0¢/0x. Equation (13) is the same as (3). Thus (10) is 
an auto-Backlund transformation for the Fisher equation. We have from 
(7), (9) and (10) that 


v=2( Ft gg) note (14) 
where ¢ satisfies (12), (13) with (11). 


Remark. The technique described above is called the Painlevé test in liter- 
ature. 


Problem 6. Consider Fisher’s equation 


du Ou 
at = Aa? t u(1 u) 
Show that 


is a traveling wave solution of this equation. Is the the solution an element 
of L2(R?) for a fixed t? 


Solution 6. 


Problem 7. Consider the nonlinear partial differential equation 
O6\* Po  (80\" Pb 0006 Pd _ 
dy) O22 \ dx) Ay? Ox Oy Oxdy 


(i) Show that this equation has an implicit solution for ¢ of the form 


«= X(¢,y,t) = f(e, thy + A(¢,t) (2) 


where f and A are arbitrary differentiable functions of ¢@ and t. 
(ii) Show that (1) may also be solved by means of a Legendre transfor- 
mation. 


(1) 


Solution 7. 


Problem 8. The sine-Gordon equation is the equation of motion for a 
theory of a single, dimensionless scalar field u, in one space and one time 
dimension, whose dynamics is determined by the Lagrangian density 


4 
L= (ut — u2) + cos (2. LL. (1) 
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Here c is a limiting velocity while m, , and py are real parameters. u,; and 
Uz are the partial derivatives of ¢@ with respect to t and a, respectively. In 
the terminology of quantum field theory, m is the mass associated with the 
normal modes of the linearized theory, while /m? is a dimensionless, cou- 
pling constant that measures the strength of the interaction between these 
normal modes. In classical theory m is proportional to the characteristic 
frequency of these normal modes. 

(i) Let 


t 
=a to —, u— muvXr (2) 
m 


L—> — 
m 


and set c= 1. Show that then the Lagrangian density becomes 


4 
m 2 


£= F(t — uz) + 2e08u) — (3) 


with the corresponding Hamiltonian density being given by 


m2 


H = sy (ui + ur — 2cosu) + py. (4) 


(ii) By choosing 


m4 


fe: (5) 


show that the minimum energy of the theory is made zero and (4) can be 


written as i 
m7 2 2 


H= Dy (uy + us, + 2(1 — cos u)). (6) 


Solution 8. 


Problem 9. Find the solution to the nonlinear partial differential equa- 


tion 
du\*\ Gu Ou Ou O?u du\*\ Gu 
(: (5) On? | “Be Of at (+ (5) ae) 
which satisfies the initial conditions (Cauchy problem) 


Ou 
ot 


(t = 0) = W(2). (2) 


Equation (1) can describe processes which develop in time, since it is of the 
hyperbolic type if 1+ (0u/dx)? — (Ou/dt)? > 0. Show that the hyperbolic 
condition for (1) implies for the initial conditions that 


1+ .a’(x) — b?(x) > 0. (3) 
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Solution 9. We simplify (1) by introducing the new variables a, 3 
t=a(a,8), t=ta,8), 2=ula(a,8),t(a,6)) = z(a,8). (4) 
Thus we seek a solution of (1) in parametric form 
r=r(a, () (5) 


where r(a, 3) = (t(a, 8), x(a, B), z(a, B)) is a vector with components ¢, x, 
and z. If we denote the scalar product of the vectors ry and rg by rjYro, 


YyYrg:= tte — @1%2 — 21292 (6) 


then (1) can be written in the following form 


r2,Dgg — 2varpDag + rZDac = 0 (7) 
where 
oe Or Sie Or eh Or (8) 
cae 8 6B’ “8 BadB 
and 
lik Lik Zik 
Digs | taeda Zee || (9) 
tg XB 2p 


Next we construct the hyperbolic solution of (7) which satisfies the initial 
conditions. The hyperbolic nature of (7) implies that 


(rarg)” — r2r3 > 0. (10) 
Equation (7) has the following equations for the characteristics 
Ga) SO: tg) = 0. (11) 


These characteristic equations together with (7) are a system of three equa- 
tions for the three functions t(a, 3), x(a, 8), and z(a,Z). It follows from 
(7), (8) and (9) that 

Dag = 0. (12) 


Equation (12) describes a linear dependence between rows of the determi- 
nant Da,g i.e., 


To, = A(a@, 8)ra + Bla, B)rg. (13) 


Taking the scalar product ra.grq and take into account that (9) and (11) 
are valid for all a, @ we have 


10 
Yo,ela =VorgpB = 3987 =0. (14) 
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With rarg 4 0 we conclude from (14) that A = B = 0. Thus we obtain 


r, = 0; rz = 0, Tog =0 


or written in components 
ar)’ (ae\” (a2) _, 
Oa Oa Oa) 
at)’ ae)" _(az\" _, 
OB Op Ogos -* 
O7t 0 Ora Or z 


aa08” Baap” Jaap” 
Thus the general solution of (15) (or (16)) is 


r(a, 8) = 11 (a) + r2(8). 


Problem 10. Consider the Korteweg-de Vries equation 


Ou 6 Ou Pu 
Ot "Oe ' Ox3 


From (1) we can derive the iteration scheme 


(+1) = a 
U = iD ; =O, 2 ces 
where 
(i) du 
qos 
# Ox 
(i) Let 


u) (x, t) = In(a — ct). 


Show that (4) converges within two steps to an exact solution 
Cc -2 
u(x,t) = =a + 2(x — ct) 


of the Korteweg-de Vries equation (1). 
(ii) Show that 
u(x,t) = (x — ct)? 
also converges to the solution (5). 
(iii) Let 
u) (x, t) = cos(a(a — ct))~*. 


(15) 


(16) 
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Show that within two steps of the iteration we arrive at 


u(x,t) = aa a 7 + 2a” tan? (a(x — ct)). (8) 


(iv) Show that by demanding that u — 0 for |x| — oo, provides a = i,/c/2 
and (5) becomes the well known soliton solution 


u(z,t) = —sech” (Se = ct)) (9) 


Solution 10. Blender J. Phys. A: Math. Gen. 24 1991 L509 


Problem 11. Solve the initial value problem of the partial differential 
equation 


with u(0, y) = y? applying a Taylor series expansion. 
Solution 11. 
Problem 12. The Zakharov-Kuznetsov equation for ion acoustic waves 


and solitons propagating along a very strong external and uniform magnetic 
field is, for a two-component plasma 


Ou Ou 0 

a + Mag + py Aw = 0 (1) 
O? oO? fond 

Acs dx2 | Oy2 | O22” (2) 


Here u is the normalized deviation of the ion density from the average. 
Exact solitonlike solutions exist in one, two and three space dimensions. 
They depend on the independent variables through the combination 


x — ct, p:=((a@—ct)? +y?)¥/?7, r:=((x— ct)? +y? + 27)¥?, (3) 
respectively, where 


Au-(c-=)u=0 (4) 


and 
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for the three cases. (i) Show that for the one-dimensional case we obtain 
the soliton solution 


u(x,t) = 3csech?(c!/? (a — ct — x9) /2). (6) 


(ii) Show that the flat soliton (6) is unstable with respect to nonaligned 
perturbations. 


Solution 12. 


Problem 13. Show that the nonlinear Schrodinger equation, 


Ow Hw 5 
im + gga t 2elul w=0, ga=+1 (1) 
has one-zone solutions 
w(x,t) = V f(O(a,t))explip(a,t)), p=pth(O), (2) 
O(a, t) := kx — wt, y= Kx — Ot (3) 


where f(@) and h(6) are elliptic functions. 


Solution 13. 


Problem 14. Consider the partial differential equation 


au)? | (Ou? _ 
Ox Oy aa 


Show that 


is a solution. 
Solution 14. 


Problem 15. We consider the diffusion equation with nonlinear quadratic 
recombination 
Ou O7u O7u Oru 


— 2 | | | 
ee Pra : Py oy : Pe ap (1) 


where D,z, Dy and D, are constants. This equation is relevant for the 
evolution of plasmas or of charge carries in solids, of generating functions. 
Find the condition on c;, c2 such that 

6 C1 


t)= 
Wea gs2 32) 5a v?/Dz + y?/Dy + 27/Dz + ca(t + to) 
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24 x? /D, + y?/Dy + 27/Dz 
a (x?/D, + y?/Dy + 27/Dz + ca(t — to))? 


(2) 
satisfies (1). 
Solution 15. 


Problem 16. The Burgers equation determines the motion of a pressure- 
less fluid subjected to dissipation 


Ou Ou Ou 


at On Ox? @) 
(i) Show that any solution v of the diffusion equation 
dv =u 
ae ae (2) 
Ot Ox 
yields a solution of the Burgers equation via the Hopf-Cole transformation 
2 Ov 
7 v Ox (3) 


(ii) Solving (1) with respect to wu we can introduce the iteration formula 


os 
yGtD = lu) —u), — ¥ =0,1,2,... (4) 
Uz 


where u; = Ou/dt etc.. Show that for u =v" we find the sequence 


yh) = (n— ee! (5) 
v Ox 
2 Ov 
Qe = 3S i 
aE fixed point (6) 


Remark. Thus the Hopf-Cole transformation is an attractor (fixed point) 
of the iteration (4). 


Problem 17. Consider the nonlinear diffusion equation 


Ou Oo Ou 
ee Ee gps | 1 
at Ox (ux) (1) 
(i) Show that this equation admits a solution of the form 
1 a \? x 
Pah tidy | ores | sleet 
u(x,t) = - oN ais (2) 


0 fort >0,|4] >1 
Ty 
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where 
v1 := (6t)'/8 (3) 


Remark. The solution exhibits a wave-like behaviour although it is not a 
wave of constant shape. The leading edge of this wave, that is where u = 0, 
is at « = x, and the speed of propagation is proportional to t7?/%. 


Solution 17. 


Problem 18. Show that the nonlinear partial differential equation (Fisher’s 


equation) : 
0 0 ee 
(Ga-lutens () 


admits the travelling wave solution 


Solution 18. 


Problem 19. Show that the coupled system of partial differential equa- 
tions 4 9 7 9 2 
EB E n n 
i—+ 2 7 =nE = E|? 1 
toe oak OR Bae oe? i”) 


admits the solution 


E(a,t) = Eo(x — st) exp (Fx ti (» =) 7 (2a) 


2r? 
n(a,t) = 5 (2b) 
cosh* \(a — st — x0) 
where 
2(1 — s)? 
Cee eee ea (2c) 


~ cosh(A(a — st — 20)) 


The soliton solution represents a moving one-dimensional plasma density 
well which Langmuir oscillations are locked. 
Solution 19. 


Problem 20. The Navier-Stokes equation is given by 


1 
— =-u-Vu+vV7u— —Vp 
p 
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where u denotes the solenoidal 


39 


V-uz=divu = 


flow velocity field, v and p are the constant kinematic viscosity and fluid 


density. Here p is the fluid pressure. We have 


1 Geb ba + 
u-Vu:= ] uy Uae ua oye + U3 a 

(e] (é) 0 
U1 Se + U2 Gy + us Ze 

Find the time evolution of 

ug _ Ou 

Oy Oz 

v:=curlu=Vxu= oan ous 

Oug _ Our 

Ou Oy 


Solution 20. Using that divu = 0 we obtain 


Ov 
—=v-Vu—-u:Vv+rV2v 

ot 

where 
a te) fa) 

1 Ge + v2 Oy 03 Be 
vais [88 + Be + 09 t 
oe die + Bt 


and analogously for u- Vv. Note that 
{v,u} :=v-Vu-—u-Vv 


maps an ordered pair of solenoidal vector fields into a solenoidal vector 
field. 


Problem 21. Show that 
Ou Ou 
eal —— 1 
ee! cons (1) 
can be linearized into 
ot as (2) 
of 
with 
H(a,t)=t, 2(x,t)=ulz,t), a(@(2,t),A2,))=2. (3) 
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Solution 21. Applying the chain rule yields 
du OUdE _ ONDE _DUDu Ou 


dt dz dt’ OFOt OF OL OF” (4) 
and _ 

du _ dude | dUdt _ Dudu _, (5) 

Ox OFOx OtOr OFOr ~~ 


Solving (4) and (5) with respect to Ou/Ot and Ou/Oz, respectively, and 
inserting into (1) we obtain (2), where we used (3). 


Problem 22. Consider the partial differential equation 


Ou Ou 
BE Bat" @) 


Show that (1) linearizes under the transformation 


ae,t)= Pe FORO. .GEGDIC ues: 


co U(s, t) 
(2) 
It is assumed that wu and all its spatial derivatives vanish at —oo. 
Solution 22. From (2) it follows that 
at _y gg.“ du OudT  OUdDt Ou (3) 
at” Cetus CCH «OH DEE! 
and 
OF : 1 Ou(s,t). _. SOU SUE) oo \ SOU 
Ot a u2(s,t) Ot ae a: Os? ae Ox “) 
and 7 
Ou OudT _ OUAt _ ~Ou (5) 
dx OX Ox OL Ox Ox 
Therefore ae ‘ 
iu UL 
a ar (6) 
Analogously 
Ou udT | WUdt _ Ou (7) 
Ot Oxrdt Odtdt at 
and hence 


OU Ou Ou Ou Ou (8) 
OrOx Ox Ot OF’ 
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From (5) we obtain 


0 (dU\  (du\* Pu (9) 
ac \dz)  \da “Oa? 
It follows that 
ude  Pudt\ (du ot Pu (10) 
6@2 Ox | Ont Oz)  \ dx roe 


Since O0t/Ox = 0 and 0z/0x = 1/u we obtain 


au du\? Ou 
aga =u (=) | wa (11) 


Inserting (1c), (5), (7) and (12) into (2) gives the linear diffusion equation 


OE” OR e) 
Problem 23. Show that the equation 
Ou 2 O74 
pee yr ees 1 
dt Aa? ) 
is transformed under the Cole-Hopf transformation 
= 10 _ 
a(E(2x,t),f(2,t))=-—, Ha2,t)=t, F2,th=x (2) 
u Ox 
into Burgers equation 
Ou Of». OU 
dt az (« | os) (8) 


Solution 23. 


Problem 24. Show that the the equation 


is linearized by the transformation 


RI 


&(a,t) = [ u(s,t)~'ds, t(a,t) =t, u(z(a, t), E(x, t)) = u(z,t). 


—oCo 
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Solution 24. 


Problem 25. Show that the equation 


du Bu 3 /du\? 3 ui 3 Au 
= + —U + —U 

Ot Ox 2\ 0x 2 02? 4° Ox 
can be derived from the linear equation 


d¢ _ Oo 
ot Ot 


and the transformation 
1 x 
o(a,t) = exp Gi u(s. ds) : 


Solution 25. 


Problem 26. Using the transformation 


aay ‘s ds 1GH2e SEED 2e aes 


—oCo u(s, t) 


the equation 
Ou 3 0°u 
sae yc asa 
Ot Ox? 
can be recast into the differential equation 


di Fu 4194/1 (du\* Pa\ _, 
dF 06 “wor\2\az) “ae ) ~ 


(ii) Show that using the Cole-Hopf transformation 


1 dv(&,t) 
OG, E)!- 30m 


u(z,t) = 


(4) 


(3) can further be reduced to the modified Korteweg-de Vries equation in 


v. Hint. We have 


OL i 1 Ou(s,f), Pe Pulethic . “Oru a 
I I. st) Ot ds = / u(s,t) ds =u 


Solution 26. 


Os Ox? 2 
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Problem 27. Show that the transformation 


a(z,)=2, eet) =ule,t), He,t) =# (1) 
Ox 
transforms the nonlinear heat equation 
Ou Ou 
a = AM 58 (2) 
into 3 3 A 
ae ee vias 
dE AR (4 IF) ) 


Solution 27. 


Problem 28. Consider the nonlinear partial differential equation 


Ou 


Au— f(u)(Vu)? +a(x, t)Vu + b(x, t) ia 


0 (1) 
where V is the gradient operator in the variables 71,...,%,, A := VV, 


f(u) and 6(x,t) are given functions, and a(x,t) is a given n-dimensional 
vector. Show that the transformation 


i (exp (- [ f(2)dz) ) ds — u(x,t) = 0. @) 


reduces (1) to the linear partial differential equation 


Av + a(x, t)Vu + 0(x, oo = 0. (3) 
Hint. From (2) we find 
“ = exp (- | - sete) a (4a) 
a = exp (- ‘i aes sete] (4b) 
- = exp (- , a Ht] ste (- | _ He) f(u) as 


Solution 28. 
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Problem 29. Consider the two-dimensional sine-Gordon equation 


Oru Oru Ou 
dx2 " Oy2 ~~ At? 


= sin(w) 
It may be regarded as describing solitary waves in a two-dimensional Joseph- 
son junction. The Lamb substitution is given by 

u(x,y,t) = 4tan7'(M(2,y,t)). (1) 


Find the equation for M. 


Solution 29. Since 
4w(1 — w? 
sin(4arctanw) = oe 


we obtain 


OM OM 07M aM \? aM \? aM \? 
2 \ on 2 
ara) (5 dy? =) au (4) (FZ) (=) aria) 


(3) 
Problem 30. We consider the nonlinear d’Alembert equation 
u = F(u) (1) 
where u = u(x), X = (%0,21,---,2n), 
oa Oo? oe? 
= erry 5 (2) 
Oxy = OXY Ox? 


and Fu) is an arbitrary differentiable function. 
(i) Consider the transformation 


u(x) = (w(x) (3) 


where w(x) and ®(w) are new unknown functions. Show that (1) takes the 
form 


d® dao 
mg BM = 
qo” + Ty? wee F(®) (4) 


wt =(55) -(o,) -~-()- 


(ii) Show that (4) is equivalent to the following equation 


d® P Zo 2b dP 
2 4 H " )-—F(®) =0 (6 
dw ( © iB) ae (FS ae) Be O18) 


where 
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where P,,(w) is an arbitrary polynomial of degree n in w, and \ = —1,0,1. 
Moreover P,, = dP,,/dw. 
(iii) Assume that ® satisfies 


PO  d®&P, 
Show that (6) takes the form 
d® P, PO 
te | be = 
oF ( lw | + (wyw" — r) =0. (8) 


(iv) Show that a solution of the system 


w = A= wpwt =r (9) 


is also a solution of (8), and in this way we obtain a solution of (1) provided 
® satisfies (7). 


Solution 30. 


Problem 31. Consider the nonlinear wave equation in one space dimen- 
sion 
Oui du 


jm em u = 0. (1) 


1 
L£(ue, Ur, t) = 5 (ue — uz) — 9(u) (2) 
with g(u) the potential function 


Qu? +1). (3) 


(ii) Definition. A conservation law associated to (1) is an expression of 


the form 
OT (u(x, t))  OX(u(a,t)) 


ot Ox 
where T is the conserved density and X the conserved flux. T and X are 
functionals of u and its derivatives. 
(ii) Show that the quantities 


1 du\? du\? ae 2 Ou Ou 
= = eee 
pS ((#) (=) gee Te a a 


=0 (4) 
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are conserved densities of (1), provided that u? —1, Ou/Ot, and Ou/Ox tend 


to zero sufficiently fast as |x| — +00. 


Hint. The corresponding fluxes are 


Ou Ou 1 du\? du\? ilar 9 
oo Nakane 7 7a Aaa ((#) a (=) rice Te. 


(6) 


T, and T, correspond to the energy and momentum densities, respectively. 


The Euler-Lagrange equation is given by 
0 (OL rs 0 (OL OL. 0 
Ot \ Our Ox \ Oux Ou 


where L is the Lagrange density (2). 


Solution 31. (i) Inserting (2) into (7) we find (1). 

(ii) Inserting T,; and Xj into (4) yields 
Ou 07u _ Ou Oru 2 du udu du Pu _ 
Ot Ot? * Ox Oxdt Ot Ox? Ot Ox dxdt 


Ou (07u Pu 3 mf 
a \ore daz * 4 


(u 


where we used (1). 


Problem 32. Consider the nonlinear partial differential equation 
Ou Ou 
Ot? = Ox? 

(i) Show that (1) admits the solution 


utuw=0. 


ux(a — vt) = tanh((7/V2)(a — vt — 2o)) = —ug(a — vt) 


where 
a= (l—v")? 


(3) 


and xo is a constant. These solutions are of kink (K) and antikink (K) 
type traveling at constant velocity v. These solutions do not tend to zero at 


infinity, but they do connect two minimum states of the potential $(u?—1)?. 


2 


(ii) The energy and momentum densities T; and T for the kink and antikink 


are obtained by substituting (2) into (5) of problem 1. Show that 


Ty (x,t) = cosh ((7/V2)(a — vt — 20)) 


(4) 


Nonlinear Partial Differential Equations 47 


To(x,t) = 1 cosh™((/V3)(e — vt — 20)). (5) 


(iii) Show that by integrating over x we find the energy and the momentum 


+00 +00 
4 ¥ f Ayu 
Bi: T,dz = -—=, PaaS Todxz = ~—~. 6 
f. mn 8/2 bg Tee o 
(iv) Show that the densities T; and T> are localized in space, in contrast 
with ux and ug. 
(v) We associate the mass 


M? := E? — P? (7) 


to (1) and the energy center 


Show that the kink and antikink solutions of (6) take the values 


8 


Me = 5: X.=vt+ 2. (9) 


Solution 32. 


Problem 33. We say that a partial differential equation described by 
the field u(a, y) is hodograph invariant if it does not change its form by the 
hodograph transformations 


z(x,y) = u(z,y), y(z,y) =%; u(z(z,y), y(z, y)) = x. (1) 


Show that the Monge-Ampere equation for the surface u(x, y) with a con- 
stant total curvature I 


07u 07 u O7u ao i Ou a du\? 2) 
Ox? Oy? Oxdy) ‘\aa)] * \dy 


is hodograph invariant. 


Solution 33. We have to show that 


Pa Oi van on wea 
Ox? OY? OXZOY 


aS 
Q|& 
ats 
Sacer 

bo 
aN 
Q| & 
ei] ei 
ees 

bo 
We 
oe 
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Applying the chain rule we have 
Ou OUOE DUO _ 


Oe Oe Oe ) Op Oe i) 
Thus oe 
U OU 

ey Ab 

Ox Ox (40) 


Here we used 
oz du Oe _ du OY_ OD 
Ox Ox’ Oy Oy’ Ox’ Oy 


Analogously ee es ae 
u OUudz udoy 


dy db Oy | Ody. se) 
Thus an 
wu Ou 
DE Ox (7b) 
For the second order deriatives we find 
Ou OU 
and Ale 
u 
a b 
ie (8b) 


It follows that 
Ou Ou du Ou OU 
Ox Ox Ox Ox Ox 


Problem 34. Show that the nonlinear equation of the Born-Infeld type 
for a scalar field u is obtained by varying the Lagrangian 


bus — uy)? (1) 


L (tes tly; es U4) =1 (1 t u t us t 


and has the following form 
ic du? _ (ou @ _ (Ou A du\* fu Pu | Ou Ou 
‘\dar) ' \dy) ° \dz ot Ox? § Oy? ° Az2 OF? 


du\? u du\? d2u du\? du du\? d2u 9 Oru du 
Ox} Ox? Oy] Oy? Oz) Oz? ot) Ot? OxOy OxOy 
Ou Ou O?u Ou Ou 0?u ’ Ou Ou O7u Ou Ou Pu Ou Ou 07 u 


Ox Ox OxOz ~ Oy Az AYOz '~ Ot Ox Bion” Ot Oy AtOy |” At Az AtOz ais 
(2) 
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Solution 34. The Euler-Lagrange equation is given by 


0 OL | a OL a 
Or Ou, | Oy Ouy | Oz0uz,  AtOu, Ou 


=0. (3) 


Problem 35. The sine-Gordon equation is the equation of motion for a 
theory of a single, dimensionless scalar field u, in one space and one time 
dimension, whose dynamics is determined by the Lagrangian density 


4 
L(ut, Ux, U) = (ul cur) 4 de cos (4) L. (1) 


Here c is a limiting velocity while m, A, and y are real parameters. u; and 
Uz are the partial derivatives of u with respect to t and x, respectively. In 
the terminology of quantum field theory, m is the mass associated with the 
normal modes of the linearized theory, while \/m? is a dimensionless, cou- 
pling constant that measures the strength of the interaction between these 
normal modes. In classical theory m is proportional to the characteristic 
frequency of these normal modes. Let 

x t mu 

x , t : 


m m i Vr 


and set c = 1. (i) Show that then the Lagrangian density becomes 


(2) 


4 
m 2 


= x ((ué —u2)+2cosu) — p (3) 


with the corresponding Hamiltonian density being given by 


L( uz, Ux, U) 


m 
H = Dy (ur + uz — 2cosu) + p. (4) 
(ii) Show that by choosing 
p= m*/d (5) 
the minimum energy of the theory is made zero and (4) can be written as 
4 


H= DD (u? + u2 + 2(1 —cosu)). (6) 


Solution 35. 


Problem 36. Show that the following theorem holds. The conservation 
law 9 


i) 
ay  (Ou/ or, Ou/Ot, w)) + pz (Ff (Ou/dx, Ou/dt, u)) =0 (1) 
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is transformed to the reciprocally associated conservation law 


SP" (Bujar’, Ou/Ot',u)} + (F"(u/da', u/at',u)} =0 (2) 


by the reciprocal transformation 


dx' (x,t) = Tdx — Fdt, (teh) SE (3) 
ak 
T’ ’ t a 4 
ROU OR) OULU St) = areas aE a) 49) 
—F(0u/0x, Ou/Odt, u) 
F' / / = Ab 
POE ONE) = a Blea Ba) oy 
(ii) Show that 
O\_« Te x0 ee eee (5) 
Ox T' Ox" Ot Tbr OF 
Solution 36. (i) From (3) and (4) we find 
dx” = T'dr' — F'dt' = T'Tdx — (T'F + F")dt = dx (5) 
if and only if 
T' (0u/Oz', Ou/Ot', u)T(Ou/Ox, Ou/Ot, u) = 1 (6) 


and 
T'(Ou/Oz', Ou/Ot', u) F (Ou/du, 0u/dt, u) + F'(Ou/dz', Ou/Ot',u) = 0. (7) 
The result follows. 


Problem 37. Consider the coupled two-dimensional nonlinear Schrodinger 
equation 


Ou O7u Ou 


tae Baa + Vays du*uu — 2wu = 0 (1a) 
Ow Ow OF. 
B Oxr2 ry. Oy? i pora (u u) =0 (1b) 


where 3, y, 6 are arbitrary constants. Consider the following transformation 
_ x cs y = 1 
2(x,y,t) =>, uz, yt) = 7, t(a,y,t)=—-, (2a) 


opis onl @ "eee ogee ae OD) 
U L,Y, — t exp 4Gt Ayt U x L,Y, YU L,Y, o] ZY, 
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w(0, 958) = syw(@(e, yt), T(E¥st), He, yt) (2c) 


Show that uw and W satisfy the same equation with the subscripts x, y, t 
replaced by @, ¥, t. 


Solution 37. 


Problem 38. Find a solution of the partial differential equation 
Ou du \? Ou 
cree 7 Ses) pp seee peenseethy fp 
Or, ( (3+) Ox ( u) 


u(a1, 22) = f(s) = f(ai+cre), s=a14+ cx 


with the ansatz 


where c is a nonzero constant. 


Solution 38. We have 


Thus 7 
df 2 ( of s 
ds (: nee se (F pay 
Integrating 
a 
fl—c+ecf 
provides 


1 
l—c+ef = 7 (a1 + cx, +6)’ 


where b is a constant of integration. Hence 


1 
u(x1,¢2) = f(x +1+ era) = 7 = 


Problem 39. An simplified analog of the Boltzmann equation is con- 
structed as follows. It is one-dimensional and the velocities of the molecules 
are allowed to take two discrete values, +c, only. Thus the distribution 
function in the Boltzmann equation, f(x,v,t), is replaced by two functions 
us(x,t) and u_(az,t) denoting the density of particles with velocity +c or 
—c, repsectively, at point x and time t. The gas is not confined, but x varies 
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over all points of the real line (—oo, 00). It is further assumed that there 
are only two types of interaction, viz., two + particles go over into two — 
particles, and vice versa, the probalitity for both processes to occur within 
one unit of time being the same number a. Then the Boltzmann equation, 
which in absence of external forces reads 


Of (x, v, £) 


BO) ev VFlxv,t) = f PvdQly—vilo(v—vi) MUA -FAl (1) 


where f’ is the final distribution, i.e., the distribution after a collision. 
Under the assumption described above (1) translates into a system of two 
equations 


Ou Ou 
=) See tra = Oe) (2) 


where c and oa are positive constants. This model is called the Carleman 
model. The Carleman model is rather unphysical. However with its aid one 
can prove almost all those results which one would like to obtain for the 
Boltzmann equation itself — as, for instance, the existence of solutions for a 
wide class of initial conditions or a rigorous treatment of the hydrodynamic 
limit. (i) Show that as for the Boltzmann equation, the H theorem holds 
for the Caleman model: The quantity 


=f (uso,t Inu+(a, t) + u_(a, t) nu_(a,t))dx (3) 


never decreases in time. (ii) Show that there exists the following generaliza- 
tions of the H theorem. Let f be concave function, defined on the half-line 
(0,00) which is once continuously differentiable. Let 


ee | (f(y) + Flu))ax (4) 


Show that : 

—Sr(u) > 0. 5 
Thus, not only does entropy never decrease, but the same is true for all 
quasientropies. 


(iii) Show that as a consequence, all Renyi entropies 


oe (ayn f (ut (e,0 Pe Cee (6) 


never decrease. In information-theretical language, this means that all sen- 
sible measures of the lack of information are nondecreasing. In other words, 
information is lost, or chaos is approached, in the strongest possible way. 
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Solution 39. 


Problem 40. For the functions u,v : [0,00) x R — R we consider the 
Cauchy problem 


Ou OU 4. 2 Ov Ov 5.9 
Ocoee ~ 
u(0,22) = u(x), v(0,) = v0(2). (10) 
This is the Carleman model introduced above. 
(i) Define 
S:=ut+y, D:=u-v (2) 


Show that S and D satisfy the system of partial differential equations 


OS OD OD Os 
at ty cece ae oy RPI | 4 
a 2 (3) 
and the conditions u > 0, v > 0 take the form 
S>0 and §S?-D?>0. (4) 


(ii) Find explict solutions assuming that S and D are conjugate harmonic 
functions. 


Solution 40. (ii) Assume that S und D conjugate harmonic functions, 
Le., 


aS aD _ aS aD _ 


Oe Oe. eon (5) 
Let z:= a+ it and f(z):=S+iD. Then 
df trys 
ee (2) +); cER (6) 
since 
df OS OD iva Buca OS 
ee Le 5 (8 D* + 2iDS +), ag = PS: (7) 


Owing to 05/0x = 0D/0t the second equation is also satsified. Let 


g(z) := af (az), aéeR (8) 
then 
WE) _ 6? (02) = (Plaz) to =—hW@re) ta’). (9) 
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For the solution 


2F(20) c=0 
_ J 2+%f(20)20 — tf (20)z 
a aq bd (20) + idet@-*) + f(20) = id d= V/ce#0 re 
(f (20) + id)et2—20 — f(zq) + id 


of (6) we only have to study the three cases c = —1,0,1. Since the problem 
is analytic we can consider 


OER, O41 (11) 
Oe’ +1 _i1- 62e-% + i20e—* sin x 


12> Ge’? —1 1 —26e-t cosa + 62e-2t * 


We can check that for |6| < 2 —1 the conditions (4) are satisfied. Thus a 
particular solution of (1) is given by 


(12) 


__ 1sgné sinh(t — log |6|) + sin x 
2 sgn@ cosh(t — log |@| — cos x 


u(t, x) (12a) 


1 sgné sinh(t — log |6|) — sinx 


v(t, 2) = 


= : 126 
2 sgné cosh(t — log |@|) — cos x 2b) 


Problem 41. The Carleman model as an approximation of the Boltz- 
mann equation is given by the non-linear equations 


ee Us 


rae +e nS o(u? — uz); (1a) 
Ou Ou_ 


(i) Show that if uz(a,0) > 0, then ux(ax,t) > 0 for all t > 0. 
(ii) Show that 


5 f us +u_)dx = 0. (2) 


The two properties have to hold, of course, in any model that is considered 
to be of some physical relevance. 
(iii) Let 


H= ic Inu, + u_Inu_)dx (3) 
R 


Show that the analogue of Boltzmann’s H-function is given by 


dH 
—< 
dt ~ e 2) 


(iv) Show that 


Fi (u4 nus + u_ In 


Nonlinear Partial Differential Equations 55 


, Inu, —u_Inu_) <0. (5) 


u_) + oo 


u+ lu, + u—Inu_ has to be interpreted as the negative entropy density 


and c(u4 mu, —u_Inu_) a 
(v) We define 


Sp f (plug) + fu) 


s the negative entropy flux. 


(6) 


where f is a concave (or convex, respectively) function. In the special case 


f(s) =—-slns 


S is the expression for entropy i.e., —H. Show that 


if f is convex. Show that 


if f is concave. 


Solution 41. 


dS; 
—_< 

a" 
dS; 


>0 
dt 


Simple computation yields 


) _ Of (uy) Ouz | Of (u_) Ou 
ay flu) + Flu) = du, Ot Dun Ot 
__ (aflus) flu) Of dus Of(u-) du- 
7 o( Ous Ou_ ) (uw — uy) —e (54 Ox Ou- Ox ) iT) 
By integration we arrive at 
ds; __ f (Of (us) af (u-) 
Wt =o / ( ae ) wa Pus (8) 


provided that u 
then OfOu is decreasing (or 
for the physical solutions (u4 


pies) Of(u_) 
Ou, Ou 


\(u2 


if f is concave and < 0, if f 


dt 


, > O0asx—d 


ar 4 


too. If f is concave (or convex, respectively), 
increasing, respectively). We therefore obtain 
eae 


(f" (us) — flu) (ue 
is convex. And consequently, 


if f is concave 
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dsr 

dt < 0, 
In the Carleman model not entropy alone increases with time, but also 
all functionals Sy (with f being concave) that could be designated “quasi- 
entropies”. 


if f is convex. 


Problem 42. The Broadwell model can then be written as 


Of, Of4 1 
Fe eR -AL) (12) 
0 
Pe A- (18) 
Of-  Of- 
fe _ OF TR ff). (16) 


The parameter € can qualitatively be understood as the mean free path. 
This system of equations serves as a model for the Boltzmann equation. 
The limit « — 0 corresponds to a vanishing mean free parth and the fluid 
regime, while € — oo approaches free molecular flow. The locally conserved 
spatial densities 


pi=fet2fotf-, pur=fy--— fr (2) 
corresponding to mass and x momentum. Show that these are governed by 
the local conservation laws 

Op  O O(pu)  O 


5 tT 5g lll) =o poo ape So (3) 


Solution 42. 


Problem 43. Consider the hyperbolic system of conservation laws 
— + —f(u) =0 (1) 


where f : R™ — R™ is continuously differentiable. A convex function 7(u) 
is called an entropy for (1) with entropy flux gq(u) if 


o nu) + Za(u) =0 () 


x 


holds identically for any smooth vector field u(x,t) which staisfies (1). 
(i) Show that (2) follows from (1) if 


v7 On OF 8G iy 
—, Ou; Our Our’ ca 
j=l 
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(ii) Show that for m = 1, every convex function n(n) is an entropy for (1) 
with entropy flux 


aw =f "af (w)dflw). (4) 


Solution 43. 


Problem 44. Consider the so-called sine-Hilbert equation 


ich (F) = —sinu (1) 
where the integral operator H is defined by 
1 oe t 
Hu(a,t) := =P f Uys) ay, (2) 
kT fico YE 


This is the so-called Hilbert transform. P denotes the Cauchy principal 
value. Let f be a continuous function, except at the singularity c. Then 
the Cauchy principal value is defined by 


pf f(x) := lim (/ f(a)da + f(e)de) 4 (3) 
=66 e7> —0o cte 
The Cauchy principal value can be found by applying the residue theorem. 
Let 
oe 
u(x,t) :=iln 4 
(ot) mtn TE " 
where 
N 
f(x,t) = [[(@-2;(¢) (5a) 
j=l 


Sx;(t) > 0 gr dy Dt NG In #Um forn~#m. (5b) 


Here x;(t) are complex functions of t and * denotes complex conjugation. 
(i) Show that 


H (F) =F mss) (6) 


follows from (2), (4) and (5). 
(ii) Show that (1) is transformed into the form 


ta eT eH) Tae =2)0) 
rer ae rrr: a (Teen Il 


gel i=t 
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(iii) Show that (7) is equivalent to the equation 


a fo ees 2) _ ay £2 
SU f) = sll? - £?) = 9). (8) 


(iv) Show that by multiplying both sides of (7) by x — x,(t) and then 
putting x = zp, we obtain 


dan 1  T[per(en(t) — 24 (6) 


id - WN MH Ay Dh ag IN: (9) 
dt 21 []ja1G4ny (tn (4) — «;(t)) 
Solution 44. 
Problem 45. Consider the Fitzhugh-Nagumo equation 
Ou Oru 
Sp = ae tu = wu a) (1) 


where a is a constant. Without loss of generality we can set —1 <a < 1. 
Insert the ansatz 


u(x,t) = f(x, t)w(z(a, t)) + g(@, t) (2) 


into (1) and require that w(z) satisfies an ordinary differential equation. 
This is the so-called direct method and z is the so-called reduced variable. 


Solution 45. Substituting (2) in (1) we obtain 


(155) — (ee. f IF) Put (413g)? )u?+ 


Ox? } dz? Ox Ox Ox? Ot) dz 


O? () O? ) 

wat n) | (4 3 99 - NU 1) aes 
(3) 

Now we must require (3) to be an ordinary differential equation for w(z). 

The procedure using the direct method is to impose that the different re- 

lationships among the coefficients of (3) to be a second order ordinary 

differential equation. However one can equally consider it acceptable to 

reduce (3) to a first order ordinary differential equation. Setting in (3) 


(24 + 1\gf —3f9" —af 4 


Oz 
g= ae (4) 
and demanding 
OP ea 102. 1/2 Oz 
S58 ye (a+1—39)5 (5a) 
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= b 
east Noa, 35,4 “dx ' da dxdt : ) 
Og Og | 
Fa? ae 1 99 Ule— 9) =9, (5c) 
(3) becomes 
dw 3 a@t+1—39,, .1/,.dw 
t +21/2 == 2) = 

dz 02z/Ox ( dz a (6) 


which is satisfied if w verifies the first order ordinary differential equation 


d 
SF + w =0 (7) 


which could be integrated at once yielding 


a 1/2 
w(2) = = (8) 


z+2z 


By combining (2a), (4) and (8) we can write the solution as 


+21/? az 
u(a,t) = Fm oe g. 


Problem 46. Consider the nonlinear Dirac equations of the form 


3 
Soy ow — Mb + Fv) =0. (1) 
op OL, 
The notation is the following 

V1 

. 4 4 _ | ye 
wiRAsCL v=] (2) 

Wa 


M is a positive constant, 
pp = (Pd, v) = vid + bade — Y3u3 — Vida (3) 


where (-,-) is the usual scalar product in C* and the y“’s are the 4 x 4 
matrices of the Pauli-Dirac representation, given by 


0) .52 I 0 k 0 Ok _ 
vos € oi Ayn = e 0 for k=1,2,3 (4) 
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where the Pauli matrices are given by 


oe ee) 


and F': R — R models the nonlinear interaction. (i) Consider the ansatz 
(standing waves, stationary states) 


(t,x) =e" u(x) (6) 


where xp = t and x = (21,272,273). Show that u: R® — C? satisfies the 
equation 


OxK 


_{ eG) 
u(x) = wir (8) 7 (8) 


Here r = |x| and (6,¢) are the angular parameters. Show that w and v 
satisfy the nonautonomous planar dynamical system 


3 

) 
iS > y*=—u- Mu+wyu+ F(au)u = 0. (7) 
k=1 


(ii) Let 


7 + = = v(F(v? — w*) — (M —w)) (9a) 
“ awlE i? =) = (i 2): (9b) 


(ii) Hint. Notice that 
tu = uyuy + U5ue2 — UZug — UZU4 (10) 


Inserting (8) yields 


du =v? — ww’. (11) 


Solution 46. 


Problem 47. Consider a one-dimensional system to describe the electron- 
beam plasma system 


One | 0 \ af 

ae OE (ne(ue + V)) =0 (1a) 
OUe OUe a e de 
Ot tue) Or Me MeNe OX (16) 
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— + (Ue + ve = 4ren;(ue + V). (1c) 


where n¢ is the density of the beam-electron fluids, m. the beam-electron 
mass, Ue the bulk fluid velocity, p = KyT.n- the particle fluid pressure, e the 
charge on an electron and € = —V@. We study the simplest case of uniform 
plasma in the absence of an external electromagnetic field, and assume 
that the positive ions are taken to form a fixed, neutralizing background 
of uniform density n; = No = const throughout the present analysis. The 
electrons move with a beam drift velocity V corresponding to the ions. 
Assume that electrostatic perturbation is sinusoidal 


1 
Ne(x,t) = No (1 + MO sin( ka) + cos(e)) ; n(t) < 3No (2a) 
0 
Ue(x,t) = u(t)(cos(kax) — sin(2kz)), E(x, t) = E(t)(cos(kx) — sin(kz)). 
(2b) 
Define 
il: ek Kole k? 
X(t — t Va A — =— met — 
(nl), YO = 5k), Zi Sew, aia TeS, p 
(3) 
Show that when (2) and (3) are substituted into (1a)-(1c), we obtain 
“ =kVX —XY +2N,Y 
dY 
=-qX +Y?—pZ 
dt q Pp 
< aN Va Vo LRG, 


Solution 47. 


Problem 48. Consider the partial differential equation 


Ou Ou Ou du dtu 
aad or a as ee =e () 


subject to periodic boundary conditions in the interval [0, Z], with initial 
conditions u(a2,0) = uo(a). We only consider solutions with zero spatial 
average. We recall that for L < 27 all initial conditions evolve into u(x,t) = 
0. We expand the solution for u in the Fourier series 


Co 


u(a,t)= S— an(t) exp(ikn2) (2) 


n=—Cco 


Are? 
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where ky, := 2na/L and the expansion coefficients satisfy 
a—y(t) = Gp(t). (3) 


Here a denotes the complex conjugate of a. Since we choose solutions with 
zero average we have ag = 0. (i) Show that inserting the series expansion 
(2) into (1) we obtain the following system for the time evolution of the 
Fourier amplitudes 


1 [oe) 
+ (kt — k2 — idk? )an + 5ikn Y= (@mOn—m + Gmansm) =0. (4) 


m=0 


(ii) Show that keeping only the first five modes we obtain the system 


d 
Set (ua — Kar + ik(Gia2 + Gay + Gya4 + Gsas)=0 (5a) 
d 
+ (2 — 815K? Jay + th(a} + 2ara3 + 2aya4 +2305) =0 (5b) 
T + (3 — 2710k* Jag + 3ik(aya2 + Gra4 + Gas) = 0 (5¢) 
d 
ae + (tua — 415° Jag + 2ik(az + 2ayag + 2a1a5) = 0 oe) 
- + (lus — 125i5k)as + 5ik(ara4 + a2a3) = 0 Po) 
where T 
k= 25 Mn s= Bi Bae (6) 


Solution 48. (i) Since 


a ye Am (t)ikmet*™* (7) 
© =—0o 
we have 
fa) co foe) ; 
us = S> > ag(t)ag(t)ikge* Mate)”, (8) 
v q=—0O p=—0o 
Furthermore 
Ou rs OR Se: & 
ge eee (9) 
and 
O7u oS 


oo yo anne (10) 
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[o.e) 
Ou 


at year (11) 
O*u — 4 ikm® 
apd = Ome. (12) 


Applying the Kronecker delta 6, to (7), (9), (10), (11), (12) and dpi¢n 
to (8) we obtain (4) where we used (3). 


Problem 49. The modified Boussinesq-Oberbeck equations are given by 
OAW el (e= | 


7 ” Ox oe) Ox Oz Oz Ox 


Ot O 
2 2 
9 R2? | 00 | 00 
Oz Ox Oz 
(1) 
where o is the Prandtl number, R the Rayleigh number, w the stream 


function, and @ a function measuring the difference between the profile of 
temperature and a profile linearly decreasing with height and time, namely 


06 OW 
Ba ae 


Oy 00 Opdd\ =e 
Or Oz Oz dx) T,; 


6=T-Tht+ a F't+ Bz. (2) 


Let L be the horizontal extenion of the convection cells, | = H/L and 
R, = r4(14+1’)3/I? the critical Rayleigh number for the onset of convection. 
One sets r = R/R.. Suppose that the temperature T} is fixed and that Tp 
increases. We define the dimensionless parameter 


7 eae eel 
Q:= —— 


T or (3) 


This quantity is a constant once the fluid and the experimental setting are 
chosen. The dimensionless temperatures appearing in the last term of the 
second of equations (1) are now expressed by 


1 : 
Beane, ee (4) 
a a 


The boundary conditions are 
w=Aw=0 atz=0, A 


6=0, atz=0, H 


Oy - 
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Show that the equation 


dX dY Z 

—_ =— =rX—-Y-XZ4 14+— 

qo eee. See aw ( +2) 
dW Z dZ 


can be derived from (1) and the boundary condition (5) via a Fourier ex- 
pansion of the function 7(z, z,t) and 6(a, z,t). 


Solution 49. The lack of the usual boundary condition stating that 
no heat flow occurs through the walls of the convection rolls, due to the 
addition of the sink term in the internal energy equation. Thus in the dou- 
ble Fourier expansion of the functions (a, z,¢) and 0(x, z,t) the unknown 
Fourier coefficients are 


v(m, n;t) = Rw(m,n;t) =: v(m, n;t) 
A(m,n;t) = RO(m, n; t) + iSO(m, n; t) =: 1(m, n; t) + i182(m,n;t). (7) 


Linear analysis of system (1) leads us towards the choice of the Lorenz like 
truncation. The Jacobian matrix of the system in the variables 1, 01, 62, 
evaluated at the origin, is block diagonal, with characteristic equation 


[80 nO? + (6 +2)? + H((o +: 1) = o(r(m,n) = 1) 


+e?(m,n)) — o((r(m,n) — 1) — (m,n))) = 0 (8) 
where we have set 
2/272 2 q(m,n)? 
q(m,n) = 1° (m7 +-n*), R-(m,n) = mp2 
R 2Rnt 
r(m,n) = Roanay Gre 7=0,1 (9) 


The simplest nonlinear system of ordinary differential equations is now 
obtained by choosing, for fixed (m,n), the four modes 


wi(m,n), O2(m,n), O1(m,n), 92(0, 2n) 
so that, after a suitable rescaling of the variables and of the time, one 


obtains (6), with b = 4n?/(m?l? + n?). 


Problem 50. Consider the complex Ginzburg-Landau equation 


Ow Pw 
Bp Ete og 


tw — (a+ icg)|w|?w. (1) 
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(i) Show that setting a = —1 in (1), and writing 
w(x,t) = R(2,t) exp(tO(a,t)) (2) 


one obtains two real equations which, after suitable linear combination and 
division by c?, can be written as 


2 2 
toate (ova S-(B))seovee) a 


1 oO 00 () 00 
R? 2 R2 =(1 2 2 1 1 (1 2 2. 
2° at ; ot eee Be " Ox AE PAE es 180) 
Here we have introduced cz := —@c; and € := 1/c,. (ii) Make an expansion 
in € of the form 
R:=Ro + PRot--, O:=6 '(0_1+¢70, +::-) (4) 
Solution 50. (i) We see that 
|w|? = R?. (5) 
Inserting (2) into (1) yields 
= (6) 
Separating out the real and imaginary part we obtain 
= (7) 


(ii) The expansion (4) becomes meaningful for sufficiently large values of 
ce, [8 = O(1)]. This first leads to the orders e~?, where we have 


0e_,\? a Je_ 
Fo ( =) = 0, >( st) =0 (8) 


respectively. Excluding Rp = 0, we obtain 00_,/0x = 0, so that O_,; only 
depends on t. Setting 


00_1 
thi= 
Oo) = = (9) 
one gets for the next orders 
oe? 
0 = Da Ro— Ro + ORR (e); (10a) 
O58 34 10 12 2 A 1 
—9, (Ro®1) = 5 Dpto (1+7)(Ro -(1-B)Ro (e). (10) 
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For 3,7 > 0, (6a) allows spatially periodic solutions 


rutet)= [pA] ((225) 200), (11) 


Here dn(u|m) is a Jacobian elliptic function that varies between (1 —m)1!/? 


and 1 with period 2K (m) [the parameter m is between 0 and 1, and K(m) 
is the complete elliptic integral of the first kind. For m — 1 the period of 
dn goes to infinity and (7) degenerates into the pulse 


(2y/6)'/*sech(y'/7r) (12) 


while for m — 0 one has small, harmonic oscillations. 


Problem 51. Consider the cubic nonlinear one-dimensional Schrodinger 
equation, 
dw Fw r 
i + + =0 1 
ia t 952 Qu|w| (1) 


where Q is a constant. (i) Show that a discretization with the periodic 
boundary conditions wj+n = w; is given by 


dw; a Wyj41 + Wi-1 — 2w; 


k 
i 7. + Qlw;|?wh =0, k=1,2,... (2) 


1 
(a) wiP wy and (bt) WY = S(wy twa). (3) 


(ii) Show that both schemes are of second-order accuracy. (iii) Show that 
in case (2a) there are first integrals, the L? norm, 


N-1 
T= >> |w;)? (4) 
j=0 
and the Hamilton function 
N-1 a Gi 
er ~ Gaee w)| 5 ul") (5) 


j=0 


The Poisson brackets are the standard ones. Thus when N = 2 the system 
is integrable. This system has been used as a model for a nonlinear dimer. 
(iv) Show that the Hamilton function of scheme (2b) is given by (h = 1) 


N-1 
. i, 4 1 * 
H=-i ) (w; (wj-1 + W741) = Q In(1 ae 7 Owiv;) (6) 
j=0 
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together with the nonstandard Poisson brackets 


{Gm Pn} = (1 P 5 QdnPn)Brnn (7) 


and 
{dm; In} = {Pm Pn} =0. (8) 


Solution 51. 


Problem 52. Consider the sine-Gordon equation in 3 + 1 dimensions 


x = sin(x) (1) 


where 


= t aye t ae AP (2) 


and x(,y,2,t) is a real valued scalar field. The sinh-Gordon equation is 
given by 


x = sinh(x) (3) 
Let 
{1 0 {01 _ {90 -i fil 
00 := 01)? O1 = 1 0)? 02: A 0 ; 03 = 0 
(4) 
be the Pauli spin matrices and 
= + 4 2 +4 2 + (5) 
Wi ag Oy OR Ob 


Let @ denote complex conjugates. Let 0 << 60 < 27,0<@< 2m, -c<A< 
co and —co < T < ow be arbitrary parameters. We set 


U := exp(i0o, exp(—idazge 7*)). (6) 
Assume that a and £ are solutions of (1) and (2), respectively. Let 
u:=a— ip. (7) 
The Backlund transformation B is then given by 
a — i8 = B(d,0,T)a (8) 


where B(@,0,7) is the Backlund transformation operator. The functions a 


and @ are related by 
1 1 
gau= sin (52) U. (9) 
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The Backlund transformation works as follows. Let a (respectively 3) be a 
solution of (1) (respectively 3)) then solve (9) for 6 (respectively a). The 
solution then solves (3) (respectively (1)). 

(i) Show for any U such that 


Cian (10) 


all solutions of (9) must be of the form (i.e. plane travelling wave with 
speed v less than one) 


ulx,y.2t)=f(n), nis ke tly +mz—wt (11) 
where k,l,m,w are real constants and 


w 


2 2 2 2 — 
ke +l +m* —w =1, ia ay ae 


(ii) Thus show that we have a Backlund transformation between the ordi- 
nary differential equations 


defined by 


Solution 52. 


Problem 53. The equation which describes small amplitude waves in a 
dispersive medium with a slight deviation from one-dimensionality is 


0 Ou Ou Ou O7u 
Da 3 t bun + st) ee og = 0. (1) 


The + refers to the two-dimensional Korteweg-de Vries equation. The — 
refers to the two-dimensional Kadomtsev-Petviashvili equation. Let (for- 
mulation of the inverse scattering transform) 


Hey eles) (2) 
Ox 


where 


K(x, 230.t) + Feast) + f K (2, s;y,t)F(s,2;y,t)ds=0 (3) 
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and F satsifies the system of linear partial differential equations 


a ik I ee a ee 
an3 | O23 Ot ”~—si“i«éi si 


=0 (4) 


where o = 1 for the two-dimensional Korteweg de Vries equation and ¢ = 1 
for the Kadomtsev Petviashvili equation. Find solutions of the form 


F(x, 2;y,t) = a(x; y, )B(z y, t) (5) 


and 
K(a,z;y,t) = L(a;y, t)G(z yt). (6) 


Solution 53. Inserting (5) and (6) into (3) yields 


a(x; y, t) 
(1+ f° a(s;y, t)3(s;y, t)ds) | (7) 


From (2) we then have the solution of (1). 


L(a;y,t) = 


Oo na 
uy) =2rgin(1+ [ alsinnslsutds) 6) 
provided functions a and ( can be found. From (4) we obtain 
da Aa OB AB 
ato atas=° ©) 
da Aa 0p 8B 
see — — —_=0. 1 
a, + a2 0, a, 72 0 (10) 


Then a and (@ admit the solution 
a(x, y,t) = exp(—la — (1?/0)y + Bt + 6) (11) 
and similarly 
B(x, y,t) = exp(—Lz + (L?/o)y + L°t + A) (12) 


where 6, A are arbitrary shifts and 1, L are constants. This form gives the 
oblique solitary wave solution of the two-dimensional Korteweg de Vries 
equation, 


u(x, y,t) = 2a?sech? (a(x + 2my — (a? + 3m?)t)) (13) 


where 


G40) m=26=m: (14) 


— 


Nl eR 


1 
2 
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For the Kadomtsev-Petviashvili equation we have o = 7, and so if we regard 
1, L as complex constants a real solution is just (13) with m — —im and 
L=L, 


Problem 54. Show that Hirota’s operators D?(f +g) and D'’(f-g) given 
in (1) and (2) can be written as 


e (-1)-)n! AI f A"-Jg 


PEL ID~ 26 Tin Ht Bet BaF" (3) 

ae = mn (1) tr-I- il n! oti f Ortm—t-J g 
Dali F)=d 2, jlm—j)! — l(n — 1)! OF Ox OtP-FAa™—F 
(4) 


Solution 54. We prove (3) by mathematical induction. The formula (4) 
can be proven in a similar way. We first try to show that 


n 2 (_1)(2-A) mn! 5 fla) A"—Jal x! 
(J -ae) Ue)-aey = Se. 6 


= gi(n— 9)! Ord Og!r-5 


 (=1)0-9) 11 8! f(e) PF g(2') ; 
= nit —s Oxi Ox!1-3 (6) 
(5) obviously holds. If we assume (5) to be true for n — 1, then 
a 9\"" Ny) = Se DOT (n= 1! 8! fw) OP Fg(2"') 
& : x) (F(@)-9@)) = De jlm—1—-j)) Oxi Oalm—1-5 
(7) 


mele ey (F(x) -9(2") 


( 
C= So (VOD (n= 1)! 7 F(x) BF g2') 


Using (7), we have 


(- a) He) -ae'y = 


Ox Ox! am; jl(n -—1— 7)! Oxi Oglh- 1-5 


CNOA (n= 1)! (O71 f(@) OIF g(2!) _ 07 F(z) O"~Ag(z') 
> ( ) 


= ji(n—1- 7)! Oxitl — Ogin—1—-3 Oxi Og !D J 
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a0 jin —1- 7)! Orit! — Oyln—1-3 
n-1 (—1)@-1!-D (n — 1)! 84 f(x) O°-Fg(z') 
=i gl(n-—1-—7)! es ee 


we (HI) 9 (n = 1)! f (x) "Fg (2') 
= 2s jg-Vi(n—Z)! Oxi Oa!r-I 


Lm (=I) (n = 1)! 0 f(w) O"~Fg(2') 
p> gn—1— 9)! O29 Oa'r-3 


3 (CDeP GD Clr PG) ote)" a2) 
=>{ j—-Wi(n—j) | oo) Oxi Oa!M—I 


IY. pasa OO) 
+ +o 


= 3 (=1)°-A) n! & f(x) O"~Fg(2') 


j= 


gin — 7)! Oxd— Oa!™m—I 


a ECC 
a eyo 


~ (=1) "nl 8 f(x) OF 9(2') 
Sd te gyb ORF. On 


+ 


Thus, (7) is true for alln = 1,2,.... Setting x’ = x we have 


n 


De(f-g) =>, 


j=0 


(<1) nl fe) Og) 
g'n— 7)! Oxi = Oars 


Problem 55. Consider the Korteweg-de Vries equation 


Ou Ou Pu 
Ot dx Ox3 


(i) Consider the dependent variable transformation 


gi Tan (eee) 
u(x,t) = a ae . (2) 


Show that f satisfies the differential equation 


OP OOF OT JOP ORY 2 
Tora Ba ot |? Oat ge +3(S4) ae 


=0. (1) 
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(ii) Show that this equation in f can be written in bilinear form 


(DzD; + Dz)(f + f) =0. 


Solution 55. 


Problem 56. The Sawada-Kotera equation is given by 


Ou Ou Ou O7u Fu Ou 
— + 45u? 1 1 
ot a Ox a ae Ox? oars om Ox? 
Show that using the dependent variable transformation 
07 In f(z, t) 
Ox? 


and integrating once with respect to x we arrive at 


u(a,t) = 2 


(DzD;i + Di)(f + f) =0. 


Solution 56. 


Problem 57. The Kadomtsev-Petviashvili equation is given by 


6 
Ox a 


O (Ou Ou | Pu Oru 
at dx * x3} ' 


Show that using the transformation 


67 In f(z, t) 
Ox? 


the Kadomtsev-Petviashvili equation takes the form 


u(a,t) = 2 


(D,D, + DE + 3D7)(f - f) =0. 


Solution 57. 


Problem 58. Consider the system of partial differential equations 


Ou Ou Ov 
O° Oa | Ox 
Ov 8 Ou Ov Ow 
a” On | Ox | Ox 
Ow Ou Ow 

+ Ou 
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where a, 3,7, c, 6 are constants. Show that the system admits the converved 
densities 


Hj =u 
Hy =v 4 5(8— che’ 

iz _(at+B u 
Ag=u+ =e -( 5) 0) 5 y#O 


Solution 58. The conserved density Ho is obvious since the first partial 
differential equation can be written as 
Ou Of1 4 - 
=~ = 7{ sau'+v). 
Ot Ox \2 


Problem 59. The polytropic gas dynamics in 1+ 1 dimensions is of the 
form 


Ou Ou |, pOp 
a on 08 
Op O 


Be Oe (pu) 


where «x is the space coordinate, t is the (minus physical) time coordinate, 
u is the velocity, p the density and [ = y — 2. Here y is the polytropic 
exponent. The constant C' can be removed by a rescaling of p. Express this 
system of partial differential equations applying the Riemann invariants 

2 


pe Pde ~ Pat, 
ERay? (x,t) # 


ry 2(z,t) = u(a,t) A 


Solution 59. We obtain (check second equation r2) 
Ory 3 1 Ory 
wn ((+a)-(+5)») & 


Org ne 3 Ore 
at = na s)n+(T+5)n) 3 


Problem 60. The Korteweg-de Vries equation is given by 


Ou 6 Ou Pu 
at On «Ox 


=0. 
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For a steady-state pulse solution we make the ansatz 


ie y= — 3 sech? (Se . ct) 


h 
where ; 
cosh(y) | 


Find the condition on b and c such that this ansatz is a solution of the 
Korteweg-de Vries equation. 


sech(y) := 


Solution 60. We find b=c. 


Problem 61. Consider the Korteweg-de Vries equation and its solution 
given in the previous problem. Show that 


[ Viwte, Sie = 


Hint. We have 


1 
[ sechi(s)as = i; cake = 2arctan(e*). 


Solution 61. 


Problem 62. Consider the Korteweg-de Vries equation 


3 
Ou Gy ot O’u 


at on Oe 


Show that 
4 cosh(2a — 8t) + cosh(4a — 64¢) + 3 
(3 cosh(a — 28t) + cosh(3a — 36t))? 


is a solution of the Korteweg-de Vries equation. This is a so-called two 
soliton solution. 


u(x,t) = 


Solution 62. 


Problem 63. Consider the one-dimensional Euler equations 


1 
Ou Ou Op 0 


ot Ox pdx | 
Op _ Op OU _ 4 
at “Oa | Pax 


Op Ou Op 
U = 


Se oS per” 
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where u(x,t) is the velocity field, p(x, t) is the density field and p(z,t) is 
the pressure field. Here t is the time, x is the space coordinate and 7¥ is the 
ratio of specific heats. Find the linearized equation around 4, p, p. 


Solution 63. We set 
u(e,t)=t%+eu'(a,t), ple,t)=Ppt+ep'(a,t), ple,t)=P+ep'(a,t). 


Inserting this ansatz into the one-dimensional Euler equations and neglect- 
ing higher order terms of O(e?) yields the linearized system 


Ou. £00". Lop. 


Be On. Bon 
Op! OP" _ ~du' 2a 
a On | Or 

/ / / 
OE sae oa UE a 


Ot a Ox Ox 


The linearized Euler equations are often used to model sound propagation. 


Problem 64. Consider the classical it Heisenberg ferromagnetic equation 


08g, 8 
Ot Ox? 
where S = (91, 92,93)", $2+$93+53 = 1 and x denotes the vector product 


The natural boundary conditions are S(x,t) — (0,0,1) as |a| — oo. 
(i) Find partial differential equation under the stereographic projection 


2u —1l+u%+v? 
a, Sa 
Q Q 
where Q=1+4+ u2 +0. 

(ii) Perform a Painlevé test. 


(iii) The Heisenberg ferromagnetic equation in the form given for u, v is 
gauge equivalent to the one-dimensional Schrodinger equation 


2 

a + 2(u? + uv?) =0 
2 

a - 2(u? + v7)u=0 


Both systems of differential equations arise as consistency conditions of a 
system of linear partial differential equations 


Op Ov 
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where ~ = (1,2)? and U and V are 2 x 2 matrices. The consistency 
condition is given by 

OU OV 

— —- — +[U,V] =0. 

a op 1 VI 


Two systems of nonlinear partial differential equations that are integrable 
if there is an invertible 2 x 2 matrix g which depends on x and ¢ such that 


ag 4 ag 
ax? ? ay 


Are the resonances of two gauge equivalent systems the same? 


U, = gqUag"+ + Y= gVog* + 


Solution 64. (i) We obtain 
Ov Pu du\? dv\? Ou Ov 
2 4 = 
aE on (() (=) ). Anon 
Ou Ov du\?* dv \* Ou Ov 
Oa om 2( (58) (=) ye Pal toe 


(ii) Inserting the ansatz 


co lo e) 
u=¢"> oud, v=o") v4 ¢ 
j=0 j=0 
into the system of partial differential equations we obtain n = m = —1. The 
resonances are r; = —1 (twice) and rg = 0 (twice). Resonances are those 


values of 7 at which it is possible to introduce arbitrary functions into the 
expansions. For each nontrivial resonance there appears a compatibility 
condition that must be satisfied if the solution is to have single-valued 
expansions. At rg = 0 we find that ug and vg can be chosen arbitrarily. 
(iii) The dominat behaviour is the same, i. n = m = —1. However, 
the resonances of the one-dimensional nonlinear Schrodinger equation are 
ry =—l, rg =0, 73 = 3, 74 = 4. 


Problem 65. The system of partial differential equations of the system 
of chiral field equations can be written as 


Ou Ou 
Ret op Le 
Ov Ov 


where u = (uj, U2, U3)2, V = (v1, v2, U3), U2 = v? = 1, J = diag(j1, jo, j3) 
is a 3 x 3 diagonal matrix and x denotes the vector product. Consider the 
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linear mapping M : R° > so(3) 


0 U3 —U2 
M(u) = —U3 0 U1 
U2 —U1 0 


where so(3) is the simple Lie algebra of the 3 x 3 skew-symmetric matrices. 
Rewrite the system of partial differential equations using M(u). 


Solution 65. First we note that 


Thus 
OM(u) , OM(u) 


ot Ox 
OM(v)  OM(v) | _ 
At re [M(v, M(Ju)] =0 


with u? = v? = 1. 
Problem 66. The Landau-Lifshitz equation 


Ou O7u 
where u = (wu, U2, U3)’, u? = 1, K = diag(k,,k2,k3) is a 3 x 3 diagonal 
matrix and x denotes the vector product. Consider the linear mapping 
M : R? = s80(3) 


0 U3 —u2 
M(u) =| -u3 0 uy 
U2 —U1 0 


where so(3) is the simple Lie algebra of the 3 x 3 skew-symmetric matrices. 
Rewrite the system using M(u). 


Solution 66. First we note that 


Thus we find 


ae + [M(u, M(u)ex] + [M(u), M(Ku)] = 0 


with u? = 1 and M(u)yz_ = 07M/02?. 
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Problem 67. Consider the nonlinear partial differential equation 


Ou _ ou 

Ot = Ox 

Let F 
u(x,t Sey 

(x, t) ra 


Find the partial differential equation for ¢. 


Solution 67. We obtain 


2 2 
20d © a ae 


Ot 5,2 — Pr 9 


Problem 68. (i) Show that the Burgers equation 


Ou _ O7u Ou 
at Oa? | Oa 


admits the Lax reprsentation 


2 


Ox 2 Ot Ox? Ox 


where u is a smooth function of x and t. 
(ii) Show that 
[T, K]b = 0 


also provides the Burgers equation, where 


o* 6) 6) O u 


Solution 68. 


Problem 69. Consider the Burgers equation 


Ot Ox? Ox” 
Consider the operator (so-called recursion operator) 


lou, u 
pees pate 
- + 3 Oe +3 
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where 3 ¢ 
=F De) =f seas. 


(i) Show that applying the recursion operator R to the right-hand side of 
the Burgers equation results in the partial differential equation 


D 


du Bu 3 Ou 320° 
dt a3 | 2" 0x2" 4" Ox" 


(ii) Show that this partial differential equation can also be derived from the 
linear partial differential equation 


dg Oe 


Ot Ax 


and the transformation 


Sah ees € [ sna) . 


Note that 


Solution 69. 


Problem 70. Find the solution of the system of partial differential equa- 
tions 


pao, Fe pono 
“8 + fq =0, “4 +9° =0. 
Solution 70. We obtain 
He.) = Te gfx.) = ————. 


Problem 71. Consider the Kortweg-de Vries equation 


Ou F Ou Ou 
Ot Nie ' Ox 


= 0. 


Show that 
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is a solution of the Kortweg-de Vries equation, where 


Ci Ck 
FE exp(—(nf + B)t= (ny +m}0)) 


ny 


A(a,t) = det (5 + 
This is the so-called N-soliton solution. 


Solution 71. 


Problem 72. Consider the partial differential equation 


Oi; Ow. Ou we 
at O22 \dr) 


Let u(x,t) = exp(u(z,t)). Find the partial differential equation for v. 


Solution 72. We obtain 


Problem 73. Consider the one-dimensional nonlinear Schrodinger equa- 
tion 


Ob 10h a 
ia +59, 7 w=0. 
Show that 
w(a,t) = 2vsech(s) exp(i¢(s, t)) 
with 


5 = v(x — C(t), o(3,t) = s+ 6(t) 


C(t) = 2ut+ Co, d(t) = 4(u? + v")t + do. 


is a solution (so-called one-soliton solution) of the nonlinear Schrérdinger 
equation. 


Solution 73. 


Problem 74. Consider the system of nonlinear partial differential equa- 
tions 
ui Ou 


a ae a => 0, 


2 2 
a “- ( 3)? 
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(i) Show that 


u(x,t) =+tanh(s/V2), v(a,t) =0 


is a solution, where s := y(a — ct). 
(ii) Show that 


u(x,t) = +tanh(/«/2s), —-v(a, t) = (1 — &)"/?sech(./«/2s) 


is a solution, where s := y(x — ct). 
(iii) Show that 


u(a,t) =+tanh(/«/2s), v(t) = —(1 — «)!/?sech(./«/28) 


is a solution, where s := y(a — ct). 


Solution 74. 


Problem 75. Consider the one-dimensional nonlinear Schrodinger equa- 
tion 


2 
oe 5S + lulu = 0, -—0 <£<0O 
Show that 
u(x,t) = 2n exp(i¢(z, t))sech(*)(z, t)) 
where 


(x,t) = —2(Ex + 2(€? — ?)t) + $0, (ax, t) = In(a + 4Et) + Yo 


is a solution of the one-dimensional nonlinear Schrodinger equation. 


Solution 75. 


Problem 76. The Landau-Lifshitz equation describing nonlinear spin 
waves in a ferromagnet is given by 


where 
S = ($1, S2,$3)", S?+$3+53=1 


and J = diag(J1, Jz, J3) is a constant 3 x 3 diagonal matrix. Show that 
Ow Ow 


— -f ake 
Ox a Ox 
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with 21 = x, x2 = —it and 
3 
i S- ZaSaVa 
a=1 
3 as 3 
. 24, 
M=i S- ZaFa Sa e7 + 2212923 S- Zy Sada 
a,B,y=1 a=1 


provide a Lax pair. Here 01, 02, a3 are the Pauli spin matrices and the spec- 
tral parameters (21, 22, 23) constitue an algebraic coordinate of an elliptic 
curve defined by 


1 
By 2B = qa — Jp), O56 = 1, 0-8: 


Solution 76. 


Problem 77. Consider the nonlinear partial differential equation 


Ou O07u Ou 
Ot Ox? 


Consider the generalized Hopf-Cole transformation 


v(t, x) = w(t)u(t, x) exp ([ dsu(s, 0) 
with 
v(t, x) 
w(t) + fj dsv(s, t) 


Find the differential equations for v and w. 


u(t, x) = 


Solution 77. We obtain 
dv @u dw _ Ov(0,t) 


Ot Ox?’ dt Ox 


Problem 78. Consider the one-dimensional nonlinear Schrodinger equa- 
tion . 9 
UL u 
i— + —— + Iclu|? =0 
“at * Ox? vee 


where x € R and c= +1. Show that it admits the solution 


u(x,t) = cape (w € (2) + ac azn) - *)) 


t t 
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Solution 78. 


Problem 79. Let a be a positive constant. The Kadomtsev Petviashvilli 
equation is given by 


d (Ou | ‘ Ou | Ou Ou 
dz Ot da | Ox3 
Consider the one-soliton solution 


u(x, y,t) = 2k?sech? (kya + koy — wt) 


where sech(z) = 1/cosh(x) = 2/(e” + e~*). Find the dispersion relation 
w(k1, kz). 


Solution 79. Inserting the one-soliton solution into the Kadomtsev Petvi- 
ashvilli equation yields 


Thus the dispersion relation is (check) 


2, ke 
w(ky, kg) = Aky + ag 
1 


Problem 80. Find the partial differential equation given by the condition 


u Ou/Oz1 = 
dct ae beeen) - 


Find a solution of the partial differential equation. 
Solution 80. We find 

=0 
This is the Monge-Ampére equation. A solution is 


Problem 81. A one-dimensional Schrédinger equation with cubic non- 
linearity is given by 
Ow h? 0? 


ah OR 9 ger, pi=ww 


where g > 0. 
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(i) Show that the partial differential equation admits the (soliton) solution 


es mv h y 2 _ my? 2u 
Veit) EoD (i h . ut)) gm cosh(y(x — vt))’ a) (1 ou 


The soliton moves with group velocity v. The phase velocity u must be 
u<vo/2. 

(ii) Show that the partial differential equation is Galileo-invariant. This 
means that any solution of partial differential equation can be mapped into 
another solution via the Galileo boost 


r—au-—Vt, w(a,t) > exp (jnv (<- 5vt)) w(t,x — Vt). 


Show that the soliton can be brought to rest. 


Solution 81. (ii) Applying the Galileo boost with V = —v provides 


Sin? h y 
s(t, =E t : 
A aes ( 2m ) gm cosh(yz) 


It solves the partial differential equation and is the soliton at rest. 


Problem 82. The KP-equations are given by 


a) (> Ou =) , S07 


dz Ot da O38) OP 
with a =i and a = —1. Show that this equation is an integrability condi- 
tion on 
OW Pw 
Iy =a— + — —0 
¥ “By = Ox? uy 


= Op (Ou * Oule',y) a = 


Solution 82. 


Problem 83. Show that the nonlinear partial differential equation 


Ou 9, O7u du\? 
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is transformated under the transformation ut u/V1 + u? into the nonlin- 
ear diffusion equation 


2 
at = 2 at u). 
x 


Solution 83. 


Problem 84. Consider the metric tensor field 
9 = 911 (41, ©2)dx1@dx14+912(41, V2)dx1@dxq+Go1(X1, L2)dx2@dx1+922(x1, L2)dxz@dxq 


with gj2 = gai and the g;, are smooth functions of 71, x2. Let det(g) = 
g11g22 — 912921 # 0 and 


>) Gil 922 912 
R= date 0g11/021 0922/0x1 Og12/0x1 
g Og11/0%2  Ogo2/Ox2 Ogi2/Ox2 


Find solutions of the partial differential equation R = 0. Find solutions of 
the partial differential equation R = 1. 


Solution 84. 


Problem 85. Consider the one-dimensional nonlinear Schrodinger equa- 
tion : 2 
Ow w 2 
Consider the ansatz 
w(a,t) = exp(twt)u(a). 


Find the ordinary differential equation for wu. 


Solution 85. We obtain 


2 
oe wu tu = 0. 


Problem 86. Consider the system of partial differential differential equa- 
tions 


2 jy 9 0 

OS gg Oi Sua), FAO Gn 0 
ot Ox? : 
0 0 ¥ 


where S = ($1, $2, $3)? and $7+S3+5? = 1. Express the partial differental 
equation using p(x,t) and p(x,t) given by 


S, = V/1—p?cos(q), So=VW1—p?sin(g), 53 =p. 
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Solution 86. We obtain 


Oq _ 1 Op 
Ot 1 — p? Ox? 
Op scO-G 
at Page 


Problem 87. Consider the two-dimensional sine-Gordon equation 
Oru Au 1 du 


One| One COR a: 


Let 

u(a1, t2,t) = 4arctan(v(21, v2,t)). 
Find the partial differential equation for v. Separate this partial differ- 
ential equation into a linear part and nonlinear part. Solve these partial 
differential equations to find solutions for the two-dimensional sine-Gordon 
equation. Note that 


sin(4a) = 4sin(a) cos(a) — 8 sin®(a) cos(a) 


sin(arctan(a)) = — cos(arctan(q@)) = = 
To Vi+a2 
and therefore 5 
sin(4arctan(v)) = See 
Furthermore 
0? —2v du \? 1 -0F 8 
Ont arctan(v) = +p (3) Tu? On?" 


Solution 87. We obtain 


Orv Atv 1 Ov dv \? dv\? 1 fav? 
2 = aye 
ery (= ? One = a (=) + (=) c (F) =e 


Thus we can separate the partial differental equation as 


Orv Ou 1 Av dv \? _ ({ ov aimee We 6318 ve 9 
T =v, T =v. 
Oxi = Ox? OF? 0x1 0X2 c? \ Ot 
The first equation is a linear wave equation. The second equation can be 
separated in two equations 


dv \? 1 (av\? dv \? 2 
=0, —]} =v. 
Ox, c2 \ Ot Ox 
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Consequently 
Ov 1 Ov Ov 
oe) = $§9U 
Ox Ve At’ Ox 2 
where s; = +1 and sg = +1 and 


v(x1,£2,t) = f(a, te”. 


This leads to Of /Ox, = (s,/c)Of/Ot with a smooth function f(x1,t) = 
f(a. — sict). As final solution we obtain 


v(a1,¥2,t) = 4arctan( f(a, — sict)e°?"?). 


Problem 88. (i) Find a non-zero vector field in R? such that 
V -curlV = 0. 
(ii) Find a non-zero vector field in R? such that 


V x curl(V) = 0. 


Solution 88. 


Problem 89. Let f:.R— R be a smooth function. Consider the partial 
differential equation 


Ou £ Of(u(x)) Ou 
Ot Ox Ox? 
Find solutions of the form u(a#,t) = ¢(a — ct) (traveling wave solutions) 


where ¢ is a smooth function. Integrate the obtained ordinary differential 
equation. 


Solution 89. We set s = «—ct. Then we obtain the ordinary differential 
equation 


dp d Bo 
a FACACLO)) aes 0 
Integrating once we arrive at 
d? 
—cb + f(6(s)) + $8 = Cy 


where C; is a constant of integration. One more integration yields 


2 5 $ 
1) = C2 +016 + $6? — F(0(s)), F(O) = i fluday 


88 Problems and Solutions 


where C2 is another constant of integration. 


Problem 90. Consider the partial differential equation (Thomas equa- 
tion) 

Oru Ou Ou  Oudu 

@udt “Ou Ot Ou Ot 
where a, b are constants. Show that the equation can be linearized with the 
transformation 


u(x,t) = —ba — at + In(v(z,t). 


Solution 90. We obtain the linear hyperbolic equation 


Problem 91. The Korteweg-de Vries equation is given by 
Ou du | Bu 


ot Oo " Oa3 0 
Setting u = O0v/Ox we obtain the equation 
Ou Ov 07u Ov a 


Oxdt Ox Ox? 


(i) Show that this equation can be dervied from the Lagrangian density 


1 D 
L= =e +(u,)? + 5 Vex)” 


where the Lagrange equation is given by 
O(OL\ 9 (A og OL OL. 0 
OV \ Ou.) On Ove) 02? Ate. J O0,.°~ 
(ii) Show the Hamiltonian density is given by 


2.4808 _ 3 1 5 
oa L=—(vz) 5 (Una). 


Solution 91. 


Problem 92. Show that the Korteweg-de Vries equation 


3 
Ou 6y oe O°u 


at 5G Os 
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admits the solution 


Show that 


[eevee = 


Solution 92. 


Problem 93. Consider the one-dimensional sine-Gordon equation 
o? oO? 
oP 7 ee + w* sin(u) = 0. 
Show that 
u+(a,t) = 4arctan(exp(d~' (a — vt — X) cosh(a)) 


u_(a,t) = 4arctan(exp(—d7'(a — vt — X) cosh(a)) 


are solutions of the one-dimensional sine-Gordon equation, where v 


ctanh(a), d= c/w. Discuss. 


Solution 93. 
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Problem 94. Let «(s,t) be the curvature and T(s,t) be the torsion with 
s and t being the arclength and time, respectively. Consider the complex 


valued function 
w(s,t) = K(s, t) exp(i ds'r(s’, t)). 
0 


Show that if the motion is described by 


0 0 0? 
Ot 7 Os Os? 


where b is the binormal unit vector, then w(s,t) satisfies the nonlinear 


Schrodinger equation 


a) oe 1 
imw+ zzw+ =|w|?w =0. 
8 2 


Ot ) 


Solution 94. 
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Problem 95. Consider the first order system of partial differential equa- 
tion 


Ou Ei Ou 
Ot Ox 
where u = (u1, U2, ug)’ and the 3 x 3 matrix is given by 


Tuy U2 U3 
U2 Uy — A 0 
U3 0 U1 — A 


where [' and A are constants. Along a characteristic curve C : 7 = x(s),t = 
t(s) for the system of partial differential equation one has 


x'(s) 
t(A — ATs) = = 
det ( AI3) = 0, r (8) 
Find the Riemann invariants. 


Solution 95. The Riemann invariants of the systems of partial differen- 
tial equations are functions R(u1,u2, ug) that are constant along the char- 
acteristic. Consequently VuR = (Ry,,Ru,,Ru;) is a left eigenvector of the 
matrix 


Vuk: (A— Als) = 0. 
From det(A — AJ3) = 0 the eigenvalue equation is 
—(A — (ur — A))(A? + (A — (+ 1)u1)A + Pui(ui — A) — r?) =0 


where r? := u3 + ug. One finds the three eigenvalues 


A=um—-A, Ay = =((P4+1)u.—A+D"?), A= 5((P+1)1—A-D"?) 


Nl re 


where D = (([ — 1)u, + A)? + 4r?.. The system of partial differential 
equations VuR-(A— Als) = 0 for the Riemann invariants is given by 


OR OR OR 
r = 
( ba “) Ou, be Ou2 i Ouz 0 
OR OR 
U2 Au, (u1 A X) Bus =0 
OR OR 
U3 Buy (uy A X) dus =0. 
With polar coordinates ug = rcos@, u3 = rsin@ we obtain 
OR OR 
T _ ——. — = 
(Tuy oun, +r yr 0 


Oo 
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For A = u; — A a solution is R = u3/ue. 


Problem 96. Consider the system of partial differential equations 


av aw _., 
ee een 
ow aw. 
Of oR ee 
DF. Bass. 2A 
ange + gVW 


Let N := V+W+4Z and J := V—W. Find ON/dt + 0J/Ox and 
OJ/Ot + OV/Ox + OW/Oz. 


Solution 96. Straightforward calculation provides 
ON sa OU as OF OV GIO 
at ° Ox Ot Ov Ox 


These are conservation laws. 


0, 0. 


Problem 97. Consider the Navier-Stokes equations 


1 
oe +(v-V)v = oy ee: 


(i) Show that in the limit v — 0 the Navier-Stokes equation are invariant 
under the scaling transformation (\ > 0) 


Pig Ni A pa 


(ii) Show that for finite v one finds invariance of the Navier-Stokes equation 
ify — Althy, 


Solution 97. 


Problem 98. (i) Show that 


k2/2 


B=w 
(cosh((ka — wt)/2))?’ 


u(x,t) = 


is a solution (solitary wave solution) of the Korteweg de Vries equation 


3 
Ou Guo OPu 


Ot dx | Ox3 =0. 


(ii) Let 


u(x,t) = ral In(F(a, t)). 
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Show that F(a, t) = 1+ e*?—#*?. 


Solution 98. 


Problem 99. Let LZ > 0. Consider the one-dimensional sine-Gordon 
equation 
Oru = d?u 
OP Ax? 
with periodic boundary conditions u(# + L,t) = u(x,t). Assume that 


+sin(u) = 0 


u(x,t) =a + u(x,t), jo(e,t)| <1 
and show that up to first order one finds 


Ou Ov 
Eo Oe 


and for u(x,t) = 6(t)e’*n® one obtains 


with k, = 27n/L. Show that all modes with 0 < k? < 1 are unstable. 
Solution 99. 


Problem 100. Study the nonlinear partial differential equation 


a @ : 
(s+) uj =—4 > Kjn exp(ur), 73S 2 


k=1 


Solution 100. 


Problem 101. Two systems of nonlinear differential equations that are 
integrable by the inverse scattering method are said to be gauge equivalent 
if the corresponding flat connections U;, Vj, 7 = 1,2, are defined in the 
same fibre bundle and obtained from each other by a A-independent gauge 
transformation, i.e. if 


Og 4 -1, OG -4 
pac V, = dV. av 1 
5) go; 1=gvog + (1) 


Uy = gU2g™* + aio 
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where g(z,t) € GL(n,R). We have 


aU, av; 7 
o. be THs (2) 
Show that ae Ay; 
2 2 7 
Be Og 1 Ua, Val = 0. (3) 


Solution 101. First we notice that 


Og! _,09 _ Og! 10g _ 
ae ane a (4) 


dz 
which follows from gg~! = I, where I is the n x n identity matrix. From 
(1) it follows that 
aU, _O a -1 HOG cal PG He OG) 4 09 a4 
=—U: tgU. —gU: 5 
Ae pees tela Ug ek ed 59 ao OO) 


OV; 0 1 =) 4 Og -1 07g 1 Og 1 Og 1 
= —V: + gV: V 6 
Be ag I a age ae a 
where we have used (3). From (4) and (5) we obtain 


OU2 -1 OV2* OU2 OV; =], 
He Oe ey! age SO ee ae)? (8) 


The commutator yields 

- Og _ a Og _ 
(U1, Vi] = [gUag7* + Bae ‘Vag ' + ard | 
= [gU2g', gV2g~"] + [QUeg7', ag "J +L) +6] 


= gU2V2g~' — gV2U2g~* = g([U2, Va])g~! 


Adding (6) and (7) results in 


Thus (3) follows. 
Both equations are integrable by the inverse scattering method. Both arise 
as consistency condition of a system of linear differential equations 


Oo Og 
Ot i U(x, t, A)¢, Or = V(a,t,A)o 
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where A is a complex parameter. The consistency conditions have the form 


OU OV 
Sets dV a: 
ae 
Let 9 9 
a -1, 9G -1 = st (OG) a 
VU, := gUag + an2 Vi := gVog° + a4 
ae aU, Vj 
1 1 
GE oe ee 


Problem 102. Consider the nonlinear Schrodinger equation in one space 


dimension ab Oe 
OW GU 26 is 1 
iE + TE + 2ly ly =0 (1) 
and the Heisenberg ferromagnet equation in one space dimension 
Os oS 
Po es baer S?=1 2 
dt Oa?” 2) 


where S = (51, 52,53)". Both equations are integrable by the inverse 
scattering method. Both arise as the consistency condition of a system of 
linear differential equations 


O® Ob 
— =U(a,t,r)® — =V(a,t,A)® 3 
ae a UlHt)®, 5° = V (x,t) (3) 
where A is a complex parameter. The consistency conditions have the form 
OU OV 
Sees oe = 4 
nae tWsVi=0 (4) 


(i) Show that On = goo. 
(ii) Show that (1) and (2) are gauge equivalent. 


Solution 102. 
Problem 103. The study of certain questions in the theory of SU(2) 


gauge fields reduced to the construction of exact solutions of the following 
nonlinear system of partial differential equations 


= 0. 


u( O7u x) dudu  Oudu | AVI, Dvd 


Duby” 0202) Oud  OF07 By Oy | bz dz 


O7u Ou Ov du  OvOu 
u — + = 2 — + -|}=0 
OyOY O2z0zZ Oy OY Oz OZ 
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au Pu Ov Ou  dvOdu 
u( ——— + — 2(——+— = 0, (1) 
OyOy O2Z0z Oy Oy O02 0z 
where uw is a real function and v and 0 are complex unknown functions of 


the real variables 7,,...,24. The quantities y and z are complex variables 
expressed in terms of 71,...,24 by the formulas 


V2y = 21 + ite, V2z := 43 — its (2) 
and the bar over letters indicates the operation of complex conjugations. 


(i) Show that a class of exact solutions of the system (1) can be constructed, 
namely solutions for the linear system 


a =) =~ - 7 =0 (3) 


where we assume that u, v, and U are functions of the variables 
ris (2yg)/? = (af + 93)? (4) 


and 23, ie., for the stationary, axially symmetric case. (ii) Show that a 
class of exact solutions of (1) can be given, where 


u=u(w), v=v(w), 5 = U(w) (5) 


where w is a solution of the Laplace equation in complex notation 


Ou O7u 
eee epee 6 
yon | O20z (6) 
and u, v and 0 satisfy 
du (du\* dv do du dv du 
=0, 2 = 0: 7 
“dw (=) dw dw “dw? dw dw (7) 


Hint. Let z= a+ ty, where x,y € R. Then 
Oe NTO oO oft 0,9 (8) 
dz 2\or dy)’ a2 2\ar dy 


Solution 103. 


Problem 104. The spherically symmetric SU(2) Yang-Mills equations 


can be written as 5 A 
Y1 2 
Ot Or On? Ae (le) 


96 Problems and Solutions 


Oy. , OP1 


OE ae cae —Aip2 + Aovyi (1b) 
0A, OA 
2 (SS) 1 hs ot) (16 


where r is the spatial radius-vector and t is the time. To find partial 
solutions of these equations, two methods can be used. The first method is 
the inverse scattering theory technique, where the [L, A]-pair is found, and 
the second method is based on Backlund transformations. 

(ii) Show that system (1) can be reduced to the classical Liouville equation, 
and its general solution can be obtained for any gauge condition. 


Solution 104. Introducing the complex function 


2(r,t) = yil(r,t) + tva(r,t) (2) 


we have from the first two equations of (1) 


Now we use the new variables 
2(r,t) = R(r, t) exp(0(r, t)). (4) 
Separating the imaginary and real parts we obtain from (3) 


_ 2 MR _ 00 damR (5) 
~ At Or” 1 Or at 


Substitution of (5) into (1c) yields 


Ao 


2 2 a 
Changing the variables 
R(r,t) =r exp(g(r, t)) (7) 
we come to the classical Liouville equation 
Lg = exp(2g) (8) 


which has the general solution in terms of two harmonic functions a(r, t) 
and b(r,t) related by the Cauchy-Riemann conditions 


exp(2g) = 4 (3) + ey (1—a?—b*)~?, (9a) 
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La = Lb=0, ar BE as (9b) 


Equation (8) was obtained with an arbitrary function (r,t). Thus (5) and 
(9) give the general solutions of problem (1) with an arbitrary function 
6(r,t) and harmonic functions a(r,t) and b(r,t), with 


R=+2r (3) + uy (a? +6? -1)7?. (10) 


Problem 105. We consider the Georgi-Glashow model with gauge group 
SU(2) broken down to U(1) by Higgs triplets. The Lagrangian of the model 


1 1 
Li — TFL FM + Dug" DY? — V(9) (1) 
where 
Fe, = 0, AG — O,A® + geare A), As (2) 
Duda = OnPa = J€ave A, be (3) 
and 
Vos (ees (4) 
4 \ 


(i) Show that the equations of motion are 
Dg Ree = — J€abe(D" bn) bc, DD" ba = (m? aa rd?) ba. (5) 


(ii) Show that the vacuum expectation value of the scalar field and Higgs 


boson mass are 


m2 


(=P = (6 


and 


My = V2XF, 


respectively. Mass of the gauge boson is M, = gF. 
(iii) Using the time-dependent t’ Hooft-Polyakov ansatz 
1— K(r,t) 


a a 1 
Ag (r, t) = 0, A; (r, t) = —~€ainTn 72 ’ alt, t) = gg 


x 
— 
3 

Ss 
YS 


where r, = Xp and r is the radial variable. Show that the equations of 
motion (5) can be written as 


yaa K=(K? +H? -1) (8a) 
" \ Or? OF a 7 
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o? o? \H? 
2 2 2:2. 
(iv) Show that with 
and Mh 
ag? (OMe 


and introducing the variables € := M,,r and r := Mt, system (8) becomes 


& OO) K(K?+H?-1) 
(aa - 3a) B= ~ 
a @ _ H(QK? + B(H? - &)) 
(saa) > ——— i 
(v) The total energy of the system F is given by 
2 
_ gb _ 
C(B) = Gia = 
i H? 1/oH H\? 1 K2H2 B 
24 T 4 2 2 24 \ 2 2\2 
[ (x K+ 5 (2) + ays Ss Be ey) a 
0 
(10) 


As time-independent version of the ansatz (3) gives the ’t Hooft-Polyakov 
monopole solution with winding number 1. Show that for finiteness of 
energy the field variables should satisfy the following conditions 


H —0, k-1 as €-0 (11) 


and 
H 6, K-00 as >. (12) 


The ’t Hooft-Polyakov monopole is more realistic than the Wu- Yang monopole; 
it is non-singular and has finite energy. 

(vi) Show that in the limit G — 0, known as the Prasad-Somerfeld limit, 
we have the static solutions, 


g 


= sinh é’ A(é) = €cothé — 1. (13) 


K(é) 


Solution 105. 


Problem 106. Consider the linear operators L and M defined by 


Lu(e,t.d) = (i + U(et,d)) west,r) 


Nonlinear Partial Differential Equations 99 


Mw(a,t,r) = (5 + Veeta)) w(a, t, A). 


Find the condition on L and M such that [L, M] = 0, where [, ] denotes the 
commutator. The potentials U(a,t,) and V(az,t, A) are typicaly chosen as 
elements of some semisimple Lie algebra. 


Solution 106. We obtain 


Problem 107. Let c > 0. Give solutions to the nonlinear partial differ- 


emtial equations 
du \? , (a : 
— —] =c. 
Or, 0x2 


Solution 107. A solution is 
u(a1, 22) = a1%1 + ave + ag 


with aj + a3 =c. 


Problem 108. ‘Two systems of nonlinear differential equations that are 
integrable by the inverse scattering method are said to be gauge equivalent 
if the corresponding flat connections U;, Vj, 7 = 1,2, are defined in the 
same fibre bundle and obtained from each other by a A-independent gauge 
transformation, i.e. if 


= Og _ - Og _ 
= 1, 09 -1 = 1, 99 -1 
U=glag +a 9, Ma=glag +9 (1) 
where g(z,t) € GL(n,R). We have 
ou, OV, 
——_ - —— U,,Vi| = 0. 2 
Ot Dx + (11, Vi] (2) 
Show that aus av. 
2 2 = 
Be Oy 1 U2 Val = 0. (3) 


Solution 108. First we notice that 


OG OG” I OG ~~ Os ay 
Oe ORY? Fes yad (4) 
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which follows from gg~! = I, where J is the n x n identity matrix. From 
(1) it follows that 
OU, O ar 


oe or 29 


OVA: 3 OG ays 4 109-1, Og 1 99 199 4 
He eo IVS Ve no he ge ao 58 (6) 


= _,09 _ O?g _, Og _,0g _ 


where we have used (3). From (4) and (5) we obtain 


OU, OV, Oy 1 AV2* _ (Uz MU\ 4 
a dc 9° a? ~%G ~9\ GE” Oe 


The commutator yields 

OG 24 
Ox! 
= [gU2g~', gV2g~"] + [gU2g7', ag") +L +6) 
= gU2V2g~* — gV2U2g~* = g([U2, Val)g~* 


Adding (6) and (7) results in 


-1 


[U1, Vi] = [gUag7* + ,gV29 


Thus (3) follows. 


Both equations are integrable by the inverse scattering method. Both arise 
as consistency condition of a system of linear differential equations 


Od Od 
Ot (x, * )d, Ox (a, 9 \o 
where A is a complex parameter. The consistency conditions have the form 
OU OV 
Oo eg a 
Let 3 3 
az -1, 99 ,-1 Pr = ee rs 
Uns gag ag Vay gvag. head 
an aU, aV, 
1 1 
C6 Ore 


Problem 109. Consider the nonlinear Schrodinger equation in one space 


dimension 
OY Ow 


7 
ta t 52 t2lyld =0 (1) 
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and the Heisenberg ferromagnet equation in one space dimension 
= SX s?=1 (2) 


where S = (5;,52,53)". Both equations are integrable by the inverse 
scattering method. Both arise as the consistency condition of a system of 
linear differential equations 
O® Og 
Ot U(a, ,A) ? Ox Via, >) ) ( ) 
where A is a complex parameter. The consistency conditions have the form 
OU OV 
pasts Ss =0 4 
apt OVI (4) 
(i) Show that $1 = g¢o. 
(ii) Show that (1) and (2) are gauge equivalent. 


Solution 109. 
Problem 110. The study of certain questions in the theory of SU(2) 


gauge fields reduced to the construction of exact solutions of the following 
nonlinear system of partial differential equations 


Ou Oru Oudu Oudu OvVOV DVD 0 
U = 0. 
OyOyY O2z0z Oy OY OzO0zZ OyOY Oz0zZ 
O?v O07 Ov Ou  OvOu 
u — + — 2 — + -|}=0 
OyOy  O0z0zZ Oy OY Oz OZ 
2a 25 = = 
ee ad ( 26, (1) 
OyOy O20z Oy Oy OzZO0z 
where uw is a real function and v and @ are complex unknown functions of 
the real variables 71,...,24. The quantities y and z are complex variables 
expressed in terms of 21,...,x4 by the formulas 
V2y := 21 + ix, V22 := 23 — ing (2) 


and the bar over letters indicates the operation of complex conjugations. 
(i) Show that a class of exact solutions of the system (1) can be constructed, 
namely solutions for the linear system 

Ov Ou Ov Ou 


6y On = eee ae (3) 
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where we assume that u, v, and U are functions of the variables 
r= (2g)? = (af + 03)? (4) 


and x3, ie., for the stationary, axially symmetric case. (ii) Show that a 
class of exact solutions of (1) can be given, where 


u=u(w), v=v(w), 0 = U(w) (5) 


where w is a solution of the Laplace equation in complex notation 


Ou O7u 
yay | O22” (6) 
and u, v and 0 satisfy 
du du\” — dv do dv. dv du 
u =, “ogee (7) 
dw? dw dw dw dw? dw dw 


Hint. Let z= a+ ty, where x,y € R. Then 
OF 2 AO Orn ti Oe (8) 
dz° 2\an dy)? a2 2\ar" “dy 


Solution 110. 


Problem 111. The spherically symmetric SU(2) Yang-Mills equations 
can be written as 


0 ) 
- = Aoy2 Ai (1a) 
0 O 
oa a = —Aiye + Aoi (1b) 
OA, OAo 
2 as 2,2 
2 (Se St) 1 ts ot) (1 


where r is the spatial radius-vector and t is the time. To find partial 
solutions of these equations, two methods can be used. The first method is 
the inverse scattering theory technique, where the [L, A]-pair is found, and 
the second method is based on Backlund transformations. 

(ii) Show that system (1) can be reduced to the classical Liouville equation, 
and its general solution can be obtained for any gauge condition. 


Solution 111. Introducing the complex function 


2(r,t) = yi(r,t) + tva(r,t) (2) 
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we have from the first two equations of (1) 


Now we use the new variables 
2(r,t) := R(r, t) exp(20(r, t)). (4) 
Separating the imaginary and real parts we obtain from (3) 


06 «~OlWMR 06 «~OlMR 
Ay ae + Or ? NS a Ot (5) 


Substitution of (5) into (1c) yields 


r’L(In R) = R? - 1, Li= sat 5p: (6) 
Changing the variables 
R(r,t) = rexp(g(r, t)) (7) 
we come to the classical Liouville equation 
Lg = exp(2g) (8) 


which has the general solution in terms of two harmonic functions a(r, t) 
and b(r,t) related by the Cauchy-Riemann conditions 


exp(2g) = 4 (#)'+ (3) (l—a? —87)-?, (9a) 


da 0b Oa Ob 
EPH cae ane 2) 


Equation (8) was obtained with an arbitrary function 0(r,t). Thus (5) and 
(9) give the general solutions of problem (1) with an arbitrary function 
6(r,t) and harmonic functions a(r,t) and b(r,t), with 


R=+2r (3+ (my) (a? +b? —1)7t. (10) 


Problem 112. We consider the Georgi-Glashow model with gauge group 
SU(2) broken down to U(1) by Higgs triplets. The Lagrangian of the model 
is 


1 a va i a a 
Lis — 7 Fi FM" + 5Dyo" DY 6" — V(¢) (1) 
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where 
Fo, = OA — O,A® + geare A), As (2) 
Duda = OnPa + J€ave Ar be (3) 
and ; 
RG 33 ee 
v(é):=-2 (¢ z =) | (4) 


(i) Show that the equations of motion are 
D,Rwve — — J€abe(D" bn) bc, DD" dba = (m? oF r$) ba. (5) 


(ii) Show that the vacuum expectation value of the scalar field and Higgs 


boson mass are 
2 


2) po me 

Me aos (6) 
and 

My = V2HF, 


respectively. Mass of the gauge boson is M, = gF. 
(iii) Using the time-dependent t’ Hooft-Polyakov ansatz 


1— K(r,t 1 
Ag (r, t) -_ 0, Ai(r, t) = —€ainTn (r ) 5 galt, t) = a ae 
Y i 


x 
— 

3 

SS 
<a 


where r, = £, and r is the radial variable. Show that the equations of 
motion (5) can be written as 


2 a a 2 2 
2 2 2 
(55 ay) ta 8 (2K? mtr a XE). (8b) 
(iv) Show that with 
gcd Mi 
PME 


and introducing the variables € := M,,r and r := M,,t, system (8) becomes 


oO? oO? K(K? + H?-1 
(Fa 7 5) = : a 7 po 
 P\.  H(2K?+B(H? -€)) 
(7 7 a) se €? a 
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(v) The total energy of the system F is given by 


= OE, 
Ne 
Tt Ey ig TSP OMS HN ody ek oy (ORE 2. Piva 
[ (weaned (F =) + yak? - 1h + + E(B 
0 
(10) 


As time-independent version of the ansatz (3) gives the ’t Hooft-Polyakov 
monopole solution with winding number 1. Show that for finiteness of 
energy the field variables should satisfy the following conditions 


H —0, k—-1 as €-0 (11) 


and 
H-€, Kk—-0 as >. (12) 


The ’t Hooft-Polyakov monopole is more realistic than the Wu- Yang monopole; 
it is non-singular and has finite energy. 

(vi) Show that in the limit G — 0, known as the Prasad-Somerfeld limit, 
we have the static solutions, 


. 


cae sinh €’ 


A(€) = €cothé — 1. (13) 


Solution 112. 


Problem 113. Consider the linear operators L and M defined by 


Lw(a,t, A) = Ge + u(e.t.a)) w(a,t, A) 


ot 


Find the condition on L and M such that [L, M] = 0, where [, ] denotes the 
commutator. The potentials U(a,t,) and V(a,t, A) are typicaly chosen as 
elements of some semisimple Lie algebra. 


Mw(a,t,r) := (‘3 + Veesta)) w(a, t, A). 


Solution 113. We obtain 
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Problem 114. Show that the Korteweg-de Vries and nonlinear Schrodinger 
equations are reductions of the self-dual Yang-Mills equations. We work on 
R* with coordinates 2% = (x, y,u,t) and metric tensor field 


g = dx ® dx — dy © dy+ du ® dt — dt © du 


of signature (2,2) and a totally skew orientation tensor €abca = €fabea]- We 
consider a Yang-Mills connection Dg := Og — Aq where the A, where the 
A, are, for the moment, elements of the Lie algebra of SL(2,C). The Ag 
are defined up to gauge transformations 


Ag: — hAgh™ — (Ogh)h™ 


where h(%_) € SL(2,C). The connection is said to be self-dual when (sum- 
mation convention) 


1 
3¢iblDe, Dal = [Da; Ds} (3) 


Solution 114. This is equivalent to the following three commutator 
equations 


[Dz + Dy, Dy] =0 
[Dz — Dy, De + Dy] + [Du, Di] =0 
[Dz — Dy, Di] =0. 


These follow from the integrability condition on the following linear system 
of equations 


Los = (D.. = Dy, + ADu)s = 0, Iys = (Dt + (Dz + Dy))s =0 


where A is an affine complex coordinate on the Riemann sphere CP! (the 
”spectral parameter”) and s is a two component column vector. We put 


Dy := Oy — A, Dy := 0, — B 
Di = OF —C, Dy, := Oy — D. 


Now require that the bundle and its connection possess two commuting 
symmetries which project to a pair of orthogonal spacetime translations 
one timelike and one null. In our coordinates these are along 0/Oy and 
0/Ou. We now restrict ourselves to gauges in which the components of the 
connection, (A, B,C, D) are independent of u and y. We also impose the 
gauge condition A+ D=0. The gauge transformations are now restricted 
to SL(2,C) valued functions of t alone under which A and B transform by 
conjugation, B —- hBh™!, etc. The equations now reduce to 


0,B=0 
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[0,, — 2A, 0; — C] =0 

20,A — [B,C] =0,B. 
These equations follow from the integrability conditions on the reduction 
of the linear system 


Los = (0, — 2A + B)s =0 


[ys = (Q—-C+0,)s = 0. 


When (2a) holds, B depends only on the variable t, so the gauge freedom 
may be used to reduce B to a normal form. The reductions are partially 
classified by the available normal forms. When B vanishes, the equations 
are trivially solveable, with the result that the connection may be put in 
the form Agdx* = A(t)d(# + y). Equation (1) is then satisfied. Thus we 
assume that B is everywhere non-vanishing. The matrix B then has just 


two normal forms 
(Ops 0 0 
1 0 


(B=«(4 21) 


We assume that the type of B is constant. In these case of type (, as t 
varies, & becomes a non-zero function of t. When B is in the Lie algebra 
of SU(2), SU(1,1) or SL(2,R), « is non-zero and is either real or pure 
imaginary. Case (a) leads to the Korteweg-de Vries equation, and case (3) 
leads to the nonlinear Schrodinger equation. 


A detailed analysis of the further reduction of the remaining equations leads 
to the following. 


The self-dual Yang-Mills equations are solved with B of type a by the 


ansatz 
2A= ( as ) 
da- —q 
—q? —2q 
( 2w —(de Zz Gye 
where 


4w = qexx — 49922 — 2q2 Ss "Ger 


a subscript x or t denotes differentiation with respect to that variable, and 


provided that q satisfies 
3 2 
PL arg 
Ot = Ox3 Ox 
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With the definition u = —q, = tr(BC) we obtain the Korteweg-de Vries 
equation 


4 + 12u 


at Ax? Ox" 

Conversely, every solution of the equations have type a, with trace(AB) 
everywhere non-zero may be reduce to this form, at worst after suitable 
co-ordinate and gauge transformations. 


Ou Ou Ou 
xr 


When tr(AB) is identically zero, the equations are explicitly solveable. 


The self-dual Yang-Mills equations are solved, with B of type @, by the 


ansatz ; i a 
= s _( ev ve. 
2a= (9 a anc = ( Ht bie 


provided w and w 
2KYt = Vex + Qua), Deby = —Won = apap 


and 2« = 1 or —i. Conversely, every solution of the equations for type @ 
may be reduced to this form. 


Chapter 3 


Lie Symmetry Methods 


Problem 1. Show that the partial differential equation 
Ou Ou 


x2 Oa 


is invariant under the transformation 


Zi (Gi, Ba €) = Lig Eo(1,%2,€) = €41 + Xp 
el - u(x1 £2) 
U(Z1(£1, £2), Z2(41, £2), €) = (1G. 


where € is a real parameter. 


Solution 1. We have 


Ou = Ou O81 | OU OZ2 = Ou, Ou 1 Gc 
dx, OF, 0%, OF, 0X, OF, OF. ~— (1 — ©)? \Azy 
Analogously 

OU OU OE, =U AT AU 

06. On On On Ors One 
Therefore 


Ou 1 Ou a )4 Ow 
0% (1—€)? \ Org a “Oar ; 
From (3) and (5) it follows that 


Ou 1 Ou Ou Ou 
OX (1 — «)? 


Dn (1 — eu) + oe e(1 — eu) 
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Inserting (2b), (3) and (6) into 


yields (1). 


Problem 2. Consider the Korteweg-de Vries equation 


Ou Ou | Ou 


Ot , ee ' Ox as 


Show that the scaling 
z=ce, tt=ct, u(a,t) =c?u'(a’,t’) 
leaves the Korteweg-de Vries equation invariant. 


Solution 2. 


Problem 3. Consider the partial differential equation 


(Ou, (Vu), u) = 0 


where ® is an analytic function, u depends on x, 21, ..., 2» and 


ae Ou Ou ae (vay? = (24), ey 


0x9 O2x4 OLn, Oxo Ox1 


du \? 
Ory) — 
Show that the equation is invariant under the Poincaré algebra 


6) fe) 0 
0x9’ 0x1’? OXn, 


Xo , Xo + x Megha 
) 


pee 
1 pay Oxo Ox 


Ox 


Solution 3. 


Problem 4. Consider the nonlinear partial differential equation (Born- 
Infeld equation) 


Ou 2 07u Ou Ou O7u Du 2 eu 
(: (zi) \e on Ot OnOt (: ($=) ae 
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Show that this equation admits the following seven Lie symmetry generators 


qe aad, gee 
44 =25, +ta Zs ue 0k 
Zo = uy, tha ; By =e +t tue. (2) 
Solution 4. 
Problem 5. Consider the Harry-Dym equation 
a _ wee =) (1) 


(i) Find the Lie symmetry vector fields. 
(ii) Compute the flow for one of the Lie symmetry vector fields. 


Solution 5. (i) The Lie symmetry vector field is given by 
0 6) 3) 

V= rly, t, a. ’ t, ay ’ t, at 2 

ne(stw) 2 bm(a,t.u)o. + 6(0,t,u) (2) 


There are eight determining equations, 


Om Om One _ Po 
Ou” Ox O} Ou =a Ou2 ” 
a6 Om _, 6 _ a _, 
Quoc On? ”—i“(itC HC CO 
ae On On om On 
3 | x 3 De pei Loge yy ales = 
aU pr OE ee ae ag ee: 


These determining equations can easily be solved explicitly. The general 
solution is 


Ne = ky thax t+ksz?, m= ko—3kyt, b= (kg+ky+2ksx)u (4) 


where ky,...,k5 are arbitrary constants. The five infinitesimal generators 


then are ‘ F 
Gi = aa? 2=% 
0 0 ee 0 
G3 =25 tun, G4 sta tun G5 = Bp 1 285, (5) 
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Thus (1) is invariant under translations (G; and G2) and scaling (G3 and 
G4). The flow corresponding to each of the infinitesimal generators can be 
obtained via simple integration. 

(ii) Let us compute the flow corresponding to G's. This requires integration 
of the first order system 


Ge... a6 dt du ~ 
Wee? lem — = 2tU (6) 


together with the initial conditions 
x(0) = 2, t(0) =t, u(0) = u (7) 


where e€ is the parameter of the transformation group. One obtains 


x = U 


*Q=qrgy HO=t W)=Goay (8) 


We therefore conclude that for any solution u = f(x,t) of (1), the trans- 
formed solution 


eee a os 
ue,8) =(0+e8)°s (2.7) (9 
will solve 
Ou 1 OP a 


Problem 6. Consider the Magneto-Hydro-Dynamics equations and carry 
out the Lie symmetry analysis. We neglect dissipative effects, and thus 
restrict the analysis to the ideal case. The equations are given by 


2 4 (v-V)ptpV-v=0 (1a) 

0 (Se +@- Vw) + V+ 5H) (H-V)H=0 (1b) 
Of + (v-V)H+HV-v—(H-V)v=0 (1e) 
V-H=0 (1d) 
g(s)swo(2)-0 


with pressure p, mass density p, coefficient of viscosity «, fluid velocity v 
and magnetic field H. 
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Solution 6. Using the first equation, we eliminate p from the last equa- 
tion, and replace it by 


a 


ap t RP(V-v) + (v- V)p =0. (2) 


If we split the vector equations in scalar equations for the vector compo- 
nents, we have a system of m = 9 equations, with p = 4 independent 
variables and g = 8 dependent variables. For convenience, we denote the 
components of the vector v by vz, vy and vz, not to be confused with partial 
derivatives of v. The variables to be eliminated are selected as follows: for 
the first 7 variables and the ninth variable we pick the partial derivatives 
with respect to t of p, vz, Vy, Uz, Hz, Hy, Hz and p. From the eighth 
equation we select OH,/Ox for elimination. We consider the case where 
k £0. We find 222 determining equations for the coefficients of the vector 
field 


a a a ,a a a 
= p-+ =~ +05 +n, 4+05-+# 
Wag he oe | oy ane Ob 
eT ee ee eee ee 
Wage Bue ag ee) One “ome e Cone 


The the components of the Lie symmetry vector fields are given by 


n” = ko os kst = key = kgz + kya 


n = kg ket kgx ky9z a kuy 


no =ka + ket + kox + kioy + kirz 


nt => ky a Kyat 
o@? =—2(kii — ki2 — ki3)p 
gp? = 2ki3p 


gp’? =ks — kgvy — kuz + (ki1 — ki2) v2 
py =ke + kgve — kiove + (Kir — kia) vy 
gp’? =kz + kove + kiovy + (k11 — ki2)vz 
o'* = ki3H, — kgHy — koH, 

oy = kygHy + kgHe — kioHe 

o'* = kigH, + koHy + kioHy. 


There is a thirteen dimensional Lie algebra spanned by the generators 


a a a a 
Gy Ot’ Go ox’ G3 Oy’ mee 

a 8 a 2a a 2a 

Oot an put Oe ag 89g a, 
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Gay Ugg + Yap, “Be He ag ~ Hope 
Co= 95 re tO "a Aya. OE 
Gio = ope be a + Hage Meg 
Cae U5; hes 25. bet tg tte 
Gio = i + 25 (vz yy in t Vz —) 
Gra = 2955 + 25 + Hegre Hy og + Heap (5) 


The Magneto-Hydro-Dynamic equations are invariant under translations 
G, through G4, Galilean boosts Gs through G7, rotations Gg through Gj, 
and dilations G; through G3. 


Problem 7. The stimulated Raman scattering equations in a symmetric 
form are given by 


Ov1 Ove Ov3 


—— = iayv9v — = 142030} — = 1430} U9. (1) 
Ox Be ns: 


The a; are real coupling constants that can be normalized to +1 and we 
have set vy = 1A}, vg = Ag, v3 = X, where Aj, Ag and X are the complex 
pump, Stokes and material excitation wave envelopes, respectively. The 
stars denote complex conjugation. Equations (1) are actually a special 
degenerate case of the full three wave resonant interaction equations. Find 
the similarity solutions using the similarity ansatz 


s(x, t) := at © similarity variable (2) 
and 


vi (a,t) = t~ = p1(s) exp(td1(s)), ve(a,t) = foe exp(—iealnt) p2(s) exp(ide(s)) 


us(¢,t) = =e" ps(s) explids(s)) (3) 


where p; > 0,0 < ¢ < 27 and € = +1. 


Solution 7. We substitute (3) into (1), introduce a phase 


b= b1 + 2 + 3 + alns (4) 
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separate out the real and imaginary parts of the equations and obtain a 
system of six real ordinary differential equations 


d 1 . d 1 
oPh ——p2p3 sin , pl Oe: — = P2P3 COs P 
ds Ss dt 8 
d 1 . d 1 
P — sp, sind, ao = —p3p1 cos@ 
ds 8 dt 8 
d dp 
“P38 —€pip2 sin ¢, p3—— = —€1p2 COs o. (5) 
ds dt 
System (5) allows two first integrals 
a 
T= pit 3, Iz = pip2p3 cos 6 — ae (6) 


These two first integrals can be used to decouple the equations. In terms 
of p; we have 


dp. 


p2 = (hh — pi)”, pssin g = —s—— (I — pt) 7? 
_.t fad a 9 I 2)—1/2 7 
p3.cos($) = mm at spi(h — 61)". (7) 


The amplitude p,; then satisfies a second order nonlinear ordinary differen- 
tial equation 


ds 


dp, dp, _ Pi dp, \" (Ia + $p7)? a 3, & 
| =-—s t{ Te ds P1i—€pi41 : 


as? ' ds I, - p? spi(1 — p?) 


Under the transformation 


(8) takes the form 


5W(W +1) 


d2W 1 1 dW\? 1dW (W-1)? B\ ¥ 
ds? a wa) (#) 8 ae 8? (aw +) +20 
(10) 


where the constants are 


12h - 212 
ans (FZ +a) : b= ae y= 2eh, 6=0. (11) 


Equation (11) is one of the six irreducible Painlevé equations, namely the 
one defining the Painlevé transcendent Py(s;a, 6, 7,6). For y = 6 = 0 
this transcendent can be expressed in terms of elementary functions. When 


W-1 
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7 #0, 6 =0, Py(s;a, 8,7,0) can be expressed in terms of the third Painleé 
transcendent Prrr. 


Problem 8. The motion of an inviscid, incompressible-ideal fluid is 
goverened by the system of equations 


— +u:Au=-—Ap (1) 


A-u=0 (2) 


first obtained by Euler. Here u = (u,v,w) are the components of the 
three-dimensional velocity field and p the pressure of the fluid at a position 
x = (2, y, Z). 

(i) Show that the Lie symmetry group of the Euler equations in three 
dimensions is generated by the vector fields 


=a tao anes, v= b5 +05. ys 
pease ee. ere 
Oz Ow Op Ot 
O 0 6) 0 O 6) O 3) 6) 
PL e oR op On On! Pe OR” On Be a Oe 
Oe er eae a ee 
Ox Oy Ou Ov Ox Oz Ou Ow 
) O O O 
te “Oy Voz Ds "Ow? i a 15 (4) 


where a,b,c, q are arbitrary functions of t. 

(ii) Show that these vector fields exponentiate to familiar one-parameter 
symmetry groups of the Euler equations. For instance, a linear combination 
of the first three fields vz + vy + V- generates the group transformations 


1 
(x,t, up) > (x+ca,t,ut ea’,p—ex-a” + 5ea-a’) 


where € is the group parameter, and a := (a,b,c). These represent changes 
to arbitrarily moving coordinate systems, and have the interesting con- 
sequence that for a fluid with no free surfaces, the only essential effect 
of changing to a moving coordinate frame is to add an extra component, 
namely, 


2 


1 W 
RET SG a:a 


to the resulting pressure. 
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(iii) Show that the group generated by vo is that of time translations, 
reflecting the time-independence of the system. (iv) Show that the next 
two groups are scaling transformations 


Si: (x, tu, p) — €X, et, u, p) 


*p). 

The vector fields rz, ry, r2 generate the orthogonal group SO(3) of simul- 
taneous rotations of space and associated velocity field; e.g., rz, is just an 
infinitesimal rotation around the x axis. The final group indicates that 
arbitrary functions of t can be added to the pressure. 


So : (x,t, u,p) > (x, et,e7u,e~ 


Solution 8. 


Problem 9. Consider the nonlinear diffusion equation 


Ou 0 1 Ou 
at ax (= =) , (1) 


Within the jet bundle formalism we consider instead of (1) the submanifold 


F(a, t, U, Ux, Ure) = Ut pas +2 —(0 (2) 


and its differential consequences 


F(a, t, U, Ug)? ++) = Ute r =0,... (3) 
with the contact forms 
6 = du — uzdx — u;dt, 0, = duy — Up, dx — Utydt,... (4) 


The non-linear partial differential equation (1) admits the Lie point sym- 
metries (infinitesimal generator) 


3) ) 
Ay = Usa Ty = ea 
3) 3) 
Sy = (—Luz — 2tur) a Vy = (tty tu). (5) 


The subscript v denotes the vertical vector fields. The non-linear partial 
differential equation (1) also admits Lie-Backlund vector fields. The first in 
the hierarchy is 


Unar VUetan 12u3 6) 
u=(% wile \z (6) 
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Find a similarity solution using a linear combination of the Lie symmetry 
vector field T,, and U. 


Solution 9. For reducing (1) we consider a linear combination of the 
vector fields T, and U, i.e. aT, + U(a € R). The equation 

(aT, + U)|6=0 (7) 
where | denotes the contraction, leads to the submanifold 


Uren  JUgUee 12u3 
t t -=0 8 
aut 43 7 ( ) 


Consequently, it follows that 
oe Une 2u2 Urn  DUgtlze 12u3 -_ 08 1 oO? 1 
es re) oe 3. fe as Be She eg | i ee i a =0 
u u u u u Ox? \ 2u Ox? \u 
(9) 
where s is a cross section s(z,t) = (x,t, u(x,t)) with js*0 = 0, js*0, = 


0,-+-. Here js is the extension of s up to infinite order. For deriving (9) 
we have taken into account the identity 


1@8u 9 dud?u 12 (), as ( 1 ) 


u® 0x3 ut OxOx? —u®_ \ Ox Ox? \ 2u? 


(10) 


and the identity given by (1). Since derivates of u with respect to t do not 
appear in (9) we are able to consider (9) as an ordinary differential equation 


of second order 2 ; aed 
=0 11 
dx3 (sa) “dx? u my 


where t plays the role of a parameter and occurs in the constant of integra- 
tion. The integration of (11) yields 


a toe Laat. (12) 
In order to determine the constants of integration C; and C2 we must 
first solve the ordinary differential equation (12), where a new constant of 
integration appears which also depends on time. Then we insert the solution 
into the partial differential equation (1) and determine the quantities C), 
Cz and C3. Equation (12) is a special case of Abel’s equation. In order to 
solve Abel’s equation we set C(t) = 2a?7/9. To simplify the calculations 
further we set C2(t) = 0. Now C3(t) is the constant of integration of the 
simplified Abel equation. Imposing the condition that u(x,t) must be a 
solution of (1) yields the following linear differential equation for C3 


—3 _ =4?C3(t) =0 (13) 
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with the solution 5 
C3(t) = kye?/94 t 


where ky is the constant of integration. Consequently, we find a similarity 
solution of the diffusion equation (1) 


3(81 — 4atkyxe(2/9)e7#}1/2 4.27 
2ax[81 — atk, xe(2/9)a7t}1/2 


u(x,t) = 


If k,; and x are positive (or both are negative), then the solution exists 
only for a finite time. If ky = 0 we arrive at the time-independent solution 
u(a,t) = K/a, where K is a constant. 


Problem 10. Consider the nonlinear partial differential equation 


Oru Ou | Ou =, 
ax3 | Vt Oe) 


where c is a constant. Show that this partial differential equation admits 
the Lie symmetry vector field 


First solve the first order autonomous system of differential equations (ini- 
tial value problem) which belongs to the vector field V, i.e. 


dt’ Fe dx’ i du! ; 
—= 3 —: =. 5 a SS" 

de de de 
From this solution of the initial value problem find the transformation be- 
tween the prime and unprime system. Then using differentiation and the 
chain rule show that the prime and unprime system have the same form. 


Solution 10. The solution of the initial value problem of the autonomous 
system of first order is given by 
t'(e) = et, x'(e) =x+ct(e® — 1), u'(€) = eu. 
Then the transformation between the prime and unprime system is 
Uia)=et 
a(t.) =a + et(e —1) 
u'(a' (a, t), t'(a, t)) = eu(z,t). 


Since Ot’ /Ot = e°, Ot’ /Ox = 0, Ox' /Ot = c(e® — 1), Ox'/Ox = 1 we have 


Ou’ — Ou! Ox' Ou! Ot! 
Ot Ox! Ot Ot Ot At 
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Thus 
Ou’ (€-1) Ow. =, Ou 
cle e=e—. 
Ox! ot! Ot 
For the first space derivative we have 
Ou! — Ou' Ox’ Ou’ Ot! Ou 


€ 


jn Oe Da Ol Oe OE” 


Thus 


— =e. 
Ox! Ox 
For the second space derivative we have 


OP OU OR ON 2 Ou Oe” 2 Oo Ot Oe Oru 


dx? Ox Ox Ox | Ox! Ox? | Ox/Ot Ox Ox *~ Ox?’ 


Analogously 
Au Ou 


= ee, 
Ox!3 Ox? 


Inserting these equations into the partial differential equation yields 


—€ Bu’ —e,/ —eE € Ou! —e, / ie Ou! 
e€ aye t ule “cle —)+ ay +e “cwe yf 
oF 3/1 / ! 
Oru of meee OMEN a pe eet 
aa +u (« c(e® — 1) 4 =) bewe “a = 0 
Thus 
OP. AcpOu. a ,Ou\ 
ae" oe 7 Bat) ~ 


Problem 11. Consider the nonlinear Schrodinger equation 


with a nonlinearity F : C — C and a complex valued w(t, 21,...,2n). We 
assume that for the given nonlinearity the energy remains bounded. We 
also assume that there exists a C!-function G : R+ — R with G(0) = 0 
such that F(z) = G’(|z|?)z for all z € C. 

(i) Show that the nonlinear Schrédinger is translation invariant. 

(ii) Show that it is phase invariant, i.e. w — e’?w. 

(iii) Show that it is Galilean invariant 


w(t,x) > eV */2e—ilvl?t/Aay (4, x—vt) 
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for any velocity v. 
(iv) Show that the mass defined by 


is a conserved quantity. 
(v) Show that the Hamilton density 


1 1 
Aw) = [ * W(t, x)? + +G(lw(t,x) Pax 
ie, 2 
is a conserved quantity. 


Solution 11. For the given problem we can work with a Hilbert space 
using the inner product 


(w,v) = de w(x)u(x) + Vw(x) - Vo(x)dx. 


Problem 12. For a barotropic fluid of index y the Navier-Stokes equation 
read, in one space dimension 


3) O —1 
rer ican (y ) Qu 
Ox 2 Ox 

Ov Ov | 2 dc Ov 

rU T Cc 

Ot Ox (y-1) Ox Ox? 
where v represents the fluid’s velocity and c the sound speed. For the class 
of solutions characterized by a vanishing presure (i.e., c = 0), the above 
system reduces to the Burgers equation. We assume for simplicity the 
value y = 3 in what follows. 


(i) Show that the velocity potential ® exists, and it is given by the following 
pair of equations 


=0 continuity equation (1a) 


=0 Euler’s equation (1b) 


Ob 0d Ov 1.4, 2 
ae aL Be 5 (v rc ). (2) 


Hint. From (2) with 


Pb 

atdx — Oxdt 
the Euler equation (1b) follows. Thus the condition of integrability of ® 
precisely coincides with the Euler equation (1b). 
(ii) Consider the similarity ansatz 


v(e,t)=f(8)F, — el#,t) = 9(s)5 (3) 
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where the similarity variable s is given by 
x 

8(x,t) := —. 4 

(x,t) a (4) 


Show that the Navier-Stokes equation yields the following system of ordi- 
nary differential equations 


f rs(7 5) SE +f -1) =0 (5a) 


2) ds 


ds? ° ds 
(iii) Show that the continuity equation (la) admits of a first integral, ex- 
pressing the law of mass conservation, namely 


ott of (2 +50 2p) +210 =f) =2%y (<4 +0). (5b) 


s’g(2f —1)=C (6) 


where C is a constant. (iv) Show that g can be eliminated, and we obtain 
a second-order ordinary differential equation for the function f. 


Solution 12. (ii) From (3) we obtain 


Ov _ dfdsx x Ov — dfdsx 1 
SE dsGeE OR” Be ede Es a) 
Oc _ dgdsx x Oc _dgdsx 1 (75) 
dt dsdtt 7’ dx dsdut | "t 
and 
uv @f (ds “cdf Psa | dfdsi | df dsl (8) 
Ox2 ds? \dx) t ' dsdx?t | dsdxt  dsdxt’ 
Since q ; d j 
S = S = —3/2 
i ge Oe (9) 
(5a) and (5b) follow. 
Problem 13. The nonlinear partial differential equation 
2 
Oru Ou ee, (1) 


021022 - Ox, 


describes the relaxation to a Maxwell distribution. The symmetry vector 
fields are given by 
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Construct a similarity ansatz from the symmetry vector field 


ae 2 ta (ng =| (3) 


Ox ee Ox1 Be 


and find the corresponding ordinary differential equation. Here c1,c2,c3 € 
R. 


Solution 13. The corresponding initial value problem of the symmetry 
vector field Z is given by 
dx’, ,  dxs du’ 


/ 
= Cy — C3% —_ = — =c3u . 4 
1> P ) de (4) 


The solution to this system provides the transformation group 


C1 Cy — C3X%1, _ ! ! 
—————e %*, x(x, u,€) = c26+%2, 2£5(x,u,€) = ue* 


(5) 
where cz # 0. Now let v2 = s/c and x; = 1 with the constant c 4 0. The 
similarity variable s follows as 


wh (x,u,e) = 2S 


(6) 


$=cry4 In 
and the similarity ansatz is 
u'(24, 2) = v(s) 
Inserting (7) into 
ul. Ou 
Ax’ Oxi, © Ox! 


leads to the ordinary differential equation 


dv dv 
coat (1 ¢) (l—v)v =0. (9) 


Problem 14. Consider the Kortweg de Vries equation 
Ou Ou | Ou 


Ot ur " x3 os 
Insert the similarity solution 
_ w(s) ee: 
MeO= Gaps § = Gals 
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where s the similarity variable and find the ordinary differential equation. 


Solution 14. We obtain 


Bw dw dw 


Problem 15. Let u be the velocity field and p the pressure. Show that 
the Navier-Stokes equations 


Ou 


ay t (a V)u— Au+ Vp =0, divu = 0 


admits the Lie symmetry vector fields 


; 0 0 ; 0 0 . 0 2 0 
Or, Ou, Or | Ou Or3 = Ous 
0 0 0 (6) 
“t Ox2 n2 Ox, = Oug ie Ou, , 
0 0 : 0 ” 0 
re 0x3 “ Ox, 2 Ouz i Oug , 
0 0 0 is 0 
Me ae "lB 2 Day 8u3 
a2 ; + + x 2 U a U 2 U 2 2 o 
Ot Oa, Oa Oa, uy un ug. Op 


Find the commutators of the vector fields and thus show that we have a 
basis of a Lie algebra. Show that \/ (a? + 73)/z3 is a similarity variable. 
Do which vector fields does it belong. 


Solution 15. 


Problem 16. The AB system in its canonical form is given by the coupled 
system of partial differential equations 


as __14\Q? 
OE = 2 On 
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where € and 7 are semi-characteristic coordinates, Q (complex valued) and 
S are the wave amplitudes satisfying the normalization condition 


aol? 
| Es er 
On 
Show that these system of partial differential equations can be written as 


a compatibility condition 


Aw12 _ A w1,2 


O€On Ono€ 


of two linear systems of partial differential equations 


O O 
ae = Faby + Cho, on eB 
“e Sis Files oo SC Ab 


Solution 16. From the compatibility condition we obtain the set of 
partial differential equation 


OF OA 
By 06 | CG — BH=0 
OG OB 
On OE 2(BF — AG) =0 
OH OC 
Dn OE 2(AH — CF) =0 
If we set 
pe ee ee 
At’ Aip On’ — Aig On 
1 1, 
F= —tp, G = 2 5 A= 79 


where yp is the spectral parameter. 


Problem 17. Show that the system of partial differential equations 


Ou, — Ou 
“Ot Oy 
Ou 10 uy, _8 Ou, 
at 3 dys | 3°) dy 
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admits the (scaling) Lie symmetry vector field 


6) 
Oug i 


0 0 0 
S= at Pa aa oe 


Solution 17. 


Problem 18. Show that the partial differential equation 


du (du? 7: O7u 
Ot \Ox) ~~ Ox? 


admits the Lie symmetry vector field 


Solution 18. 


Problem 19. Consider the Kuramoto-Sivashinsky equation 


Ou OF. O7u pow 
Ot ey "Bx? “Oat 


=0. 
Show that the equation is invariant under the Galilean transformation 


(u, x,t) (ut+c,x2—ct,t). 


Solution 19. 


Problem 20. Consider the nonlinear partial differential equation 
Pu (F =) 
+u ae =0 


where c is a constant. Show that the partial differential equation admits 
the Lie symmetry vector fields 


O 0 
“= 5p Y= op 
O O O O O 
Va = 3ta, + (x + 2ct) =, Va=ta, t+ cto” hun. 
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Solution 20. Consider the vector field V4. The autonomous system of 
first order differential equation corresponding to the vector field is 


dt’ ‘i da! 1! du! ; 
= = ae 


de’ de ’— de 


The solution of the initial value problem is 


This leads to the transformation 


t'(x,t)=et 
x (x,t) =(e§ —l1)ct+a2 
UE (Bb), 6 st) = SU, 0) 


Then since O0t'/Ox = 0, Ox’ /Ox = 1 we have 
Ou’ Ou’ Ox’ Ou’ OF 


be Bl On Oh Oe © Bx 
_ duu 
Sg. Oa 


Ou’ Ou’ Ot! — Ou’ Ox - OU 
| =e 


at ot’ at | Aa’ dt at 


_ oul. OW ec 1) Ou 
~ BH * Oa! Ce OE 
Oru OP ut On! Ou’ Ot! Ou 


On? Ba2 Oz. Bix On Ox? 
-sOPU -: 2070 
Oa! . Ox2 


Ou! OPu Os" a 
@x2 O28 Ox Ox8 
Ou! Pu 
~ O28 © Oa3" 


Inserting these expression into the partial differential equation yields 


3,7 / / / 
gee See ey (F gee 1)c4 ete ) = 
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Thus 


with e~* £0. 


Problem 21. Consider the stationary incompressible Prandtl boundary 
layer equation 

Ou Ou Pu = Ou -OFu 

On? = On OnOE —OE - ANODE 
Using the classical Lie method we obtain the similarity reduction 


u(E,n) = €°y(ax), = x = n€* + FE) 


where f is an arbitrary differentiable function of €. Find the ordinary 
differential equation for y. 


Solution 21. Differentiation and applying the chain rule provides 


3 d? d 2 
oY = py (28 1)( 4 


dx x2 dx 


which, in the special case 3 = 2, is the Chazy equation. 


Problem 22. Show that the Chazy equation 


d’y d?y dy \” 
dx3 2y dx? @ dx 
admits the vector fields 


6) ) 6) 2 O ) 

a — -yx 2ry +6 

Ox’ "Oa Voy ” Ox Qzry By 
as symmetry vector fields. Show that the first two symmetry vector fields 
can be used to reduce the Chazy equation to a first order equation. Find 
the commutators and classify the Lie algebra. 


Solution 22. 


Problem 23. Show that the Laplace equation 
Oe 0? oO” 
=0 
(a = Oy? a) 4 


admits the Lie symmetries 
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0 0 0 0 0 
Moz => — 2a EP a 


0 re) 
K,= ps et Ky = —2yD—r? 
where r? := x? + y? + 2?. 


Solution 23. 


Problem 24. Find the Lie symmetry vector fields for the Monge-Ampére 
equations 


Ot? Ox? Oxdt 
where K = +1, K =0, K =-1. 


O7u OP u ( O7u ) 7 


Solution 24. 


Problem 25. Let u = (ui, us, ug), V = (v1, V2, v3)? and uf +u3+u2 = 1, 
vi + v3 + v3 = 1. Let x be the vector product. The O(3) chiral field 
equations are the system of partial differential equations 


where R is a 3 x 3 diagonal matrix with non-negative entries. Show that 
this system of partial differential equations admits a Lax pair (as 4 x 4 
matrices) L, M, i.e. [L,M] =0. 


Solution 25. 


Problem 26. Consider the Schrédinger-Newton equation (MKSA-system) 
h2 
—S—AG + UG = pt 
m 
AU = 41Gm? ||? 


where m is the mass, G the gravitational constant and ys the energy eigen- 
values. The normalization consition is 


| (a1, 22, 23)|?daidr2dx3 = 1. 
R3 
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Write the Schrodinger-Newton equation in dimension less form. Then find 
the Lie symmetries of this system of partial differential equation. 


Solution 26. 


Problem 27. Show that the one-dimensional wave equation 


102u—s- 0?'u 
2 Ot? Ox? 


is invariant under 


(an) (Santa) cane) (4) 
u'(x'(a,t), (2, t)) = u(a,t). 


Solution 27. 


Problem 28. Study the system of nonlinear partial differential equation 


ee & ; 
(a + 5) tis = 45° Kye exp(ue) 
1 


l= 


with j = 1,2 and K is the Cartan matrix 


Solution 28. 


Problem 29. The partial differential equation 


is called the Tzitzeica equation. Applying u(r, y) = e“( show that the 
equation takes the form 
O?v Ov Ov 3 


"Oxdy Ox Oy aa 


Find the Lie symmetries of both equations. 


Solution 29. 
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Problem 30. Find the Lie symmetry vector fields for the one-dimensional 
telegraph equation 


Pu im) 


ot Ox? 
where c? = 1/(LC), a= G/C, 8 = R/L and G is the conductance. 


+ (a4 aye taBu=c 


Problem 31. Find the Lie symmetry vector fields of the Thomas equation 


0?u Ou ; gee Ou Ou 
andi | at "On ‘Ox Ot 


= 0. 


Problem 32. The Harry Dym equation is given by 
Ou 3 0°u 
pasar yc pees 
Ot 0x3 

Show that it admits the Lie symmetry vector fields 


6) ) 
Va= 5 


) ) 6) ) ) 3) 
ee oe Sag —, Ve=27?— +2ru—. 
Ber pg PAT rrp Oy ae OG 
Is the Lie algebra spanned by these generators semi-simple? 
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